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Abstract 


The  characteristics  of  several  types  of  heated  sensors  used 
for  measuring  oceanic  turbulence  have  been  examined.  The  author 
measured  the  steady  heat  flux  from  glass  coated  ellipsoidal  microbead 
thermistors,  computed  numerically  the  unsteady  heat  flux  from  an  ideal 
constant  temperature  flat  plate,  and  analyzed  the  steady  flow  calibration 
data  of  a paralene-c  coated  plate  thermistor  as  well  as  the  steady  and 
unsteady  calibration  data  of  two  conical  constant  temperature  hot  film 
anemometers.  In  order  to  understand  the  behaviour  of  probes  alone,  thermal 
models  of  these  sensors  have  been  developed.  The  models  incorporate  a 
Nusselt  number  governed  steady  heat  flux  from  the  wetted  surface  and  the 
thermal  effects  of  a surface  coating  and  a supportive  substrate.  Derived 
functional  relationships  between  the  steady  heat  flux  and  the  flow  rate 
agree  favourably  with  the  available  calibration  data.  The  quasi-steady 
sensitivity  of  these  probes  when  used  as  anemometers  or  thermometers  as 
well  as  their  signal  contamination  by  temperature  or  velocity  are 
calculated  using  the  functional  heat  flux  relationships.  The  substrate 
and  the  coating  reduce  the  sensitivity  to  temperature  and  to  speed  as 
well  as  the  ratio  of  speed-to-temperature  sensitivity. 

The  response  of  sensors  is  not  governed  by  the  Nusselt  number 
when  the  boundary  layer  is  unsteady.  The  unsteady  response  of  a sensor 
to  velocity  oscillations  is  governed  by  its  unsteady  viscous  boundary 
layer  and  may  increase  with  increasing  frequency  over  some  frequency 
ranges.  The  response  bandwidth  is  wider  for  velocity  than  for  tempera- 
ture. The  ratio  of  unsteady  temperature-to-veloc i ty  sensitivity  is 
highest  at  zero  frequency.  Frequency  response  calibrations  methods 
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1 

I . I ntroduct i on 

This  thesis  discusses  some  salient  features  of  heated  anemometry 
and  thermometry  in  water  with  emphasis  on  oceanic  applications.  The 
discussion  separates  the  hydrodynamic  character i sti cs  of  probes  from 
their  associated  electronics.  Many  pertinent  features  of  a probe's 
behaviour  can  now  be  deduced  unobscured  by  the  characteristics  of  the 
probe's  electronics.  The  three  appendices  in  this  thesis  are  manuscripts 
that  have  been  submitted  to  refereed  journals  and  deal  with  three 
different  but  related  aspects  of  the  subject  of  thermometry  and 
heated  anemometry.  This  thesis  i nter-re I ates  and  summarizes  the  contents 
of  the  three  manuscripts. 
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2.  Background 

The  direct  observation  of  turbulence  is  important  to  oceano- 
graphy. Observations  of  the  mixing  of  momentum,  salinity  and  temperature 
between  water  masses  could  clarify  the  details  of  frontal  dynamics 
(Gargett,  1978)  and  the  role  of  turbulence  in  general  circulations. 

The  measurement  of  turbulence  on  scales  small  enough  to  resolve  the 
dissipation  spectrum  is  presently  the  only  practical  way  to  estimate 
the  vertical  Reynolds  shear  stress  in  open  ocean  shear  flows  (Crawford 
and  Osborn,  1978).  The  correlation  of  orthogonal  velocity  fluctuations 
in  river  estuaries  could  yield  the  internal  stress,  an  important  but 
so  far  unmeasured  term  in  the  balance  of  forces. 

An  understanding  of  the  hydrodynamics  of  the  employed  sensors 
is  crucial  to  the  successful  interpretation  of  turbulence  data.  In- 
variably a sensor  must  be  mounted  on  a platform  such  as  a towed  body 
(Gargett,  1978;  Grant,  Stewart  and  Moilliet,  1962),  a self  Dropelled 
vehicle  (Grant,  Moilliet  and  Vogel,  1968),  a free  falling  package 
(Osborn,  1978),  a tethered  falling  package  (Mamorino  and  Caldwell,  1978) 
or  a platform  attached  to  the  bottom  or  the  shore.  The  motions  of  the 
platform  are  also  important  for  data  i nterpretation  but  are  beyond  the 
scope  of  this  thesis. 

A diagrammatic  summary  of  turbulence  sensors  that  have  been 
used  for  oceanic  measurements  is  shown  in  Figure  I.  A summary  of 
salinity  sensors  has  been  given  by  Gregg  and  Pederson  (1978).  Cylindrical 
film  and  wire  sensors,  absent  from  Figure  I,  have  so  far  only  seen 
use  in  laboratories  (Tan-atichat,  Nagib  and  Pluister,  1973  and  Gust, 

1976). 
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The  following  character i st ics  of  thermometers  and  anemometers 
used  in  water  must  be  examined  to  determine  their  suitability: 

(1)  environmental  adaptability 

(2)  electrical  insulation  from  the  fluid 

(3)  spatial  resolution 

(4)  frequency  response 

(5)  sensitivity  to  temperature  and  velocity 

Environmental  adaptability  includes  the  requirements  of  mechanical 
strength,  ridigidy  and  longevity  as  well  as  the  ability  to  shed 
contaminants  like  plankton  and  dirt.  Choosing  a sensor  that  meets  the 
requirements  of  environmental  adaptability  is  an  art  as  much  as  it  is 
a science.  For  example  the  Ocean  Mixing  Group  of  the  Institute  of  Ocean 
Sciences  at  Patricia  Bay,  Canada,  has  found  that  their  cone  shaped  hot 
film  sensors  remain  relatively  free  from  contaminants  and  can  be  success- 
fully back  flushed,  in  situ,  to  remove  larger  plankton  such  as  jelly-fish. 

The  sensing  element  of  heated  anemometers  and  thermometers 
is  a conductor  with  a temperature  coefficient  of  electrical  resistance. 

The  sensing  element  can  be  a metal -oxide  core  or  a thin  metal  film. 
Suitable  electronics  measure  the  temperature  and  the  heating  of  the 
sensing  element.  Spurious  signals  are  introduced  by  conduction  through 
sea  water  if  the  sensing  element  is  not  electrically  insulated. 

The  spatial  resolution  of  a probe  is  primarily  governed  by 
its  size.  The  turbulence  cannot  be  fully  resolved  if  its  spatial  scale 
is  comparable  to  or  smaller  than  the  probe.  A detailed  discussion  of  the 
spatial  resolution  of  probes  is  outside  the  scope  of  this  thesis.  A 
spatial  structure  of  wavelength  X advected  past  a probe  at  speed  U has, 
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for  the  probe,  an  apparent  frequency  U/x.  If  U/X  exceeds  the  frequency 
resolution  of  the  probe  then  there  is  an  apparent  spatial  response 
limitation  caused  by  frequency  response  limitations.  Frequency  response 
is  discussed  in  Chapter  4. 

Determining  sensitivity  to  temperature  or  velocity  is  a 

fundamental  problem  of  thermometry  and  heated  anemometry.  Temperature 

sensors  and  heated  anemometers  are  simultaneously  sensitive  to  both 

temperature  and  velocity.  Thermometers  have  a small  amount  of  internal 

heating  because  of  the  electric  current  passing  through  them  and 

consequently  have  some  sensitivity  to  speed.  The  level  of  internal 

heating  that  can  be  tolerated  is  not  clear.  Gregg,  Meagher,  Pederson 

and  Aagaard  (1978)  write: 

"A  third,  more  obvious,  approach  [to  increase 
the  sensitivity  of  thermistors]  would  be  to 
increase  the  current  through  the  thermistors; 
however,  by  increasing  the  power  dissipation, 
this  increases  the  ventilation  noise  [sensitivity 
to  speed].  As  neither  the  velocity  spectra  nor 
the  coefficients  for  computing  the  cooling  have 
been  we  I I known,  we  have  fol lowed  the  conserva- 
tive position  of  Cox,  et  a I (1969)  in  maintain- 
ing  a low  power  dissipation,  i.e.  4 to  13  pw 
[10~6  watts].  With  further  information  it  may 
be  possible  safely  to  increase  the  power." 

The  contents  of  the  square  brackets  have  been  added  for  clarity.  The 

heat  flux  from  grossly  heated  sensors  responds  to  variations  of  the 

ambient  fluid's  speed  as  well  as  its  temperature.  Increasing  the 

sensor's  temperature  increases  the  ratio  of  speed-to-temperature 

sensitivity.  Eventually  the  formation  of  bubbles  and  a shortening 

of  the  operating  life  limit  the  allowable  overheat.  Gargett  (1978) 

comments  on  the  problems  of  limited  overheat: 


t 
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” ...  a second  problem  arises  because  the 
hot  film  is  not,  i n 


fact,  very  hot  (measure- 
ments were  taken  with  an  overheat  of  20  C,  a 


value  chosen  as  a compromise  between  velocity 
sensitivity  and  film  corrosion,  both  of  which 
increase  with  increasing  overheat),  and  as 
such  it  is  sensitive  to  temperature  as  well 
as  velocity.  A single  test  of  the  hot  film 
probe  in  the  plume  tank  used  to  calibrate 
cold  film  temperature  probes  [Fabula,  1968] 
suggests  that  this  temperature  sensitivity 
may  be  higher  than  would  be  expected  at  the 
high-frequency  end  of  the  velocity  dissipation 
range,  high  enough  to  be  a significant  problem 
in  most  regions  of  the  main  thermocline,  where 
microscale  temperature  fluctuations  almost 
always  accompany  dissipative  velocity  events. 
Until  this  problem  receives  further  attention, 
the  velocity  data  from  our  tows  in  the  main 
thermocline  remain  suspect." 


3.  Steady  Heat  Transfer 
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Much  work  has  been  published  on  the  steady  heat  transfer 
from  ideal  bodies,  that  is,  probes  with  a constant  temperature  wetted 
surface.  Real  probes  have  a coating  that  insulates  the  sensing  element 
and  some  have  a substrate  that  conducts  heat  away  from  the  sensing 
element.  However,  the  concepts  of  heat  transfer  from  ideal  probes  can 
be  fruitfully  applied  to  real  probes  used  for  oceanic  work  to  deduce 
their  steady  and  guasi-steady  characteristics. 


3.1  Ideal  Probes 

The  rate  of  heat  transfer  from  a constant  temperature  wetted 
surface  (an  ideal  probe)  to  its  fluid  environment  depends  on  (Hinze, 

1975) 

( 1 ) the  flow  velocity 

(2)  the  temperature  difference  between  the  surface  and  the  ambient 
fluid  (over-heat  A9 ) and 

(3)  the  physical  properties  of  the  fluid. 

The  rate  of  heat  transfer  from  a given  constant  temperature  wetted  surface 
is  character i zed  by  a dimensionless  parameter  called  the  Nusselt  number 
(Nu).  The  dependence  of  the  Nusselt  number  on  other  dimensionless 
groups  of  importance  may  be  expressed  by  (Hinze,  1975) 


where 


Nu 


{%.  Pr,  Gr, 


A! 

T ' 


C 46'  Z'  0} 
P 


ql_3S  A9 
Gr  = ^ — 5 — 


( I ) 


A9=  temperature  difference  between  the  wetted  surface  and 
the  ambient  fluid  (overheat), 


I 
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L = a relevant  dimension  of  the  probe, 
k = thermal  diffusivity  of  the  fluid, 
v = kinematic  viscosity  of  the  fluid, 
p = density  of  the  fluid. 


0^  = specific  heat  of  the  fluid. 


8 = coefficient  of  thermal  expansion, 
g = gravitational  acceleration, 

U = speed  of  the  fluid, 

T = absolute  temperature  of  the  fluid, 

Z = a parameter  group  representing  the  effect  of  geometry, 

= angle  between  the  flow  direction  and  the  axis  of  the  probe 
(if  it  has  one) . 

The  term  U2/C  A9  represents  the  relative  importance  of  flow  induced 


heating  (stagnation  Temperature,  U2/CM  which  eguals  'vIMOOO  °C  for 


typical  fluid  speeds  of  100  cm/sec.  The  stagnation  temperature  is 
possibly  important  for  thermometry  (A9  n,  10" 3 °C)  but  not  for  anemometry 


where  A9  ,v  I0°C.  Collins  and  Williams  (1959)  showed  that,  for  wires  in 


air,  buoyancy  effects  (thermal  convection)  may  be  neglected  whenever 

3 


G < R 
r e 


If  the  above  condition  also  holds  for  water  and  other  shapes  then 
buoyancy  effects  are  small  whenever 


U~ 


500  — > I in  units  of  cm/sec  and  C 
A9 


a condition  that  is  always  well  satisfied  in  oceanic  work.  The  term 
A9/T  is  important  if  the  overheat  is  sufficient  to  change  either  the 
Reynolds  number  or  the  Prandtl  number.  For  water  both  the  Prandtl 
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number  and  the  kinematic  viscosity  depend  strongly  on  temperature 
(table  I),  consequently  A0/T  may  be  an  important  dimensionless  group. 
However,  the  author  is  unaware  of  any  attempt  to  relate  the  Nusselt 
number  to  A9/T  in  water.  For  anemometry  in  water  the  Nusselt  number 
can  be  expressed  as 

Nu  = f (Re,  Pr,  y2-  » Z,  *}  (2) 

and  further  discussion  is  limited  to 

Nu  = f {Rg,  Pr,  Z}  (3) 

when  Re>>l  and  I s Pr  s 10. 


Kinematic  viscosity 

T C cm2/sec  Prandtl  Number 


0 

1.787  x I0~2 

13.44 

5 

1.514 

10 

1.304 

9.45 

15 

1.138 

8.  13 

20 

1.004 

7.07 

25 

0.894 

30 

0.802 

5.49 

35 

0.725 

40 

0.659 

4.34 

50 

0.554 

60 

0.474 

3.01 

70 

0.414 

80 

0.366 

2.25 

90 

0.327 

100 

0.295 

1.78 

TABLE  I.  The  kinematic  viscosity  and  Prandtl  number  of  pure  water 
after  Batchelor  (1976). 


The  Nusselt  number  for  a sphere  and  a cylinder  is  known  from 
observations.  The  Nusselt  number  has  been  derived  analytically  and 
empirically  for  wedges,  including  a flat  plate,  and  for  cones.  The 
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/ 


Nusselt  number  of  a constant  temDerature  sphere  is  (Kramers,  1946) 

Nu  = 2.0  + 1.3  Pr0*15  + 0.66  Pr0*31  R * 

e 


(4) 


where  the  Reynolds  number  is  based  on  the  sphere's  diameter.  The  rate 
of  heat  transfer  to  the  fluid  is 


0 = nDX  f A9Nu 


(5) 


where  X^  is  the  conductivitv  of  the  fluid  at  the  sphere's  surface  and  D 
is  the  diameter  of  the  sphere.  The  Nusselt  number  for  two-dimensional 
flow  over  a wedae  ( Fa  I kner-Skan  flow)  has  been  tabulated  by  White  (1974) 
for  a wide  range  of  prandtl  numbers  and  wedge  angles.  A flat  plate  is 
a wedge  of  zero  angle  and  has  a Nusselt  number  of  0.664  Pr*" 


R z. 
e 


The  Nusselt  number  of  a two-dimensional  stagnation  area  is  .715  Pr 

I 

R 1 . The  rate  of  heat  transfer  from  both  faces  of  a wedge  is 
e 

0 = bXfA9Nu 


.57 


(6) 


where  b is  the  width  of  the  wedge  across  the  flow. 

There  is  a transformation  due  to  Mangier  (1948),  valid  for 
both  incompressible  and  compressible  laminar  flow,  which  converts  the 
axi-symmetr ic  boundary  layer  equations  of  a cone  into  the  equations  of 
an  equivalent  two-dimensional  wedge  flow.  The  heat  flux  from  any  axi- 
symmetric  body  can,  in  principle,  be  obtained  from  the  heat  flux  of 
an  equivalent  two-dimensional  wedge. 

Hinze  (1975)  gives  an  extensive  discussion  on  the  Nusselt 
number  of  a cylinder.  All  findings  are  not  completely  in  agreement 
but  a relationship,  claimed  to  be  satisfactory  for  many  gases  and 


liquids,  is  (Kramers,  1946) 


Nu  = 0.42  Pr0,20  + 0.57  Pr’/3  R 1/2 

e 


(7) 


L J 


The  heat  flux  from  a cylinder  of  length  b is 


0 = irbXfA9Nu. 


The  quasi -steady  sensitivity  of  the  heat  flux  to  speed  for 
an  ideal  probe  is 

30  _ 90  _ , A0,  , ^ 9Nu  p 

9l)  " U 9R  ( U Xf  "b  9R  Re 

e be 

where  the  factor  irb  is  used  for  cylinders,  ttD  for  spheres  and  the 
factor  b for  wedges  and  cones.  The  quasi -steady  sensitivity  to 
temperature  oscillations  of  the  ambient  fluid  is 
90  . , 770  . . 

aTSoT  = xf  ( *b  > Nu 

D 

The  ratio  of  the  quasi -steady  speed-to-temperature  sensitivity  is 

90/9U  / M.  ) My. 

9O/9(A0)  1 U ' Nu  9R 

e 

i 

This  ratio  is  A0/2U  if  Nu  « Rg2  which  is  the  case  for  wedges  and  cones 
but  is  less  than  A9/2U  for  cylinders  and  spheres.  For  very  large 
Reynolds  numbers  the  ratio  of  the  quasi-steady  speed-to-temperature 
sensitivity  for  cylinders  and  spheres  appraoches  A0/2U  asymptotically. 
Thus  the  maximum  value  of  the  ratio  of  quasi-steady  speed-to-temperature 
sensitivity  of  heated  sensors  never  exceeds  A0/2U. 


3.2  Real  Probes 

In  principle,  if  the  Nusselt  number  of  a probe  is  known  and 
if  its  temperature  and  its  heat  flux  can  be  measured  then  the  fluid's 
speed  can  be  uniquely  determined.  Two  difficulties  with  using  the 
Nusselt  number  of  an  ideal  probe  for  a real  probe  are  (aside  from 


geometric  imperfections) : 


p ip 
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(1)  All  probes  are  electrically  and  hence  thermally  insulated  from  the 
fluid.  The  wetted  surface  is  cooler  than  the  electrically  heated 
element  and  may  not  be  at  a constant  or  uniform  temperature  even 
if  the  heated  element  is  at  a constant  and  uniform  temperature, 

(2)  Cones,  wedges  and  plates  may  have  a substrate  that  is  not  completely 

i enveloped  by  the  heated  element.  If  so,  there  is  a heat  loss  into 

I 

I 

the  substrate.  Similarly  cylinders  and  spheres  may  have  a heat 
loss  into  their  supports. 

In  principle,  the  heat  flux  from  any  real  probe  can  be  calculated  by 
solving  the  coupled  temperature  equations  of  the  boundary  layer  and  the 
solid  probe  but,  in  practice,  this  would  entail  an  inordinate  amount 
of  effort.  The  heat  flux  from  any  probe  can  also  be  determined 
empirically  but  this  by  itself  would  shed  little  light  on  the  para- 
meters governing  the  heat  flux  unless  there  is  insight  to  help  interpret 
the  results. 

In  practice,  only  the  average  temperature  of  the  sensing 
element  and  the  heat  flux  from  it  are  known.  The  temperature  of  the 
wetted  surface  is  not  directly  assessable.  Thermal  models  of  real 
probes  that  incorporate  the  concept  of  a Nusselt  number  governed  heat 
flux  from  their  wetted  surfaces  and  that  also  include  the  thermal  effects 
of  the  coating  and  the  substrate  have  been  developed  (appendices  B and  C). 

The  temperature  of  the  wetted  surface  is  variable  but  is  assumed  to  be 

I 

uniform.  With  these  models  a theoretical  functional  relationship 
between  the  heat  flux  and  the  flow  can  be  derived  without  any  knowledge 
of  the  probe's  electronics.  These  relationships  take  the  form 
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— = h(R  *>  (8) 

0 e 

and  agree  favourably  wi+h  the  available  calibration  data  of  hot  film 
probes  (appendix  B)  and  heated  thermistors  (appendix  C). 

Regardless  of  qeometry,  a constant  temperature  surface  has 
a heat  flux  that  is  either  a proportional  or  a linear  function  of  U1 . 
Tan-Atichat,  et  al  (1973)  found  it  impossible  to  fit  the  heat  flux 

I 

versus  UJ  calibration  data  of  their  0.015  cm  diameter,  glass  coated, 
cyl indr i cal,  hot  film  probe  to  a linear  relationship.  They  attempted 
the  more  general  linearization 

= A + BUn  (9) 

A0 

and  claimed  a best  fit  when  A = 0 and  n = 0.23.  Equation  9 is  often 
fitted  to  calibration  data  from  hot  wires  in  air  but  A is  never  zero 
and  the  exponent  n aDnears  to  be  restricted  to  0.45  to  0.52.  There 
is  no  physical  model  to  justify  equation  9 when  n does  not  equal  j. 

The  steady  calibration  data  of  two  conical  hot  film  probes  (appendix  B) 

I 

also  show  that  the  heat  flux  is  not  a proportional  function  of  U2  , 

I 

as  expected  from  equation  6,  or  even  a linear  function  of  U2.  Similarly, 
the  calibration  data  of  a paralene-c  coated  plate-shaped  thermistor, 
a paralene-c  coated  ellipsoidal  microbead  thermistor  and  two  glass 
coated  microbead  thermistors  (appendix  C)  show  that  the  heat  flux  is 
not  a linear  function  of  U2.  The  heat  flux  relationship  does  not  take 
the  form 

~ = h(R  *)  = [A  + B R *]_1 

0 e e 

(equation  9 with  n=i)  because  of  the  thermal  effects  of  the  coating  and 
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the  substrate.  The  heat  flux  relationships  are 


= A + 


for  a coated  plate  shaped  thermistor. 


M_  - a + — 
Q A Nu 


for  a coated  ellipsoidal  microbead  thermistor  and 


— = C + 
9 


A + R 2 
e 


for  a coated  conical  film  sensor  on  a substrate.  The  oarameters, 

A,  B and  C depend  on  the  probe's  geometry,  size,  coating  thickness, 
the  thermal  conductivity  of  the  fluid,  the  coating  and  the  substrate, 
and  the  fluid's  Prandtl  number.  The  use  of  a third  parameter  for  the 
conical  probe  model  reflects  the  presence  of  the  substrate. 

With  a knowledge  of  the  function  h the  quasi-steady  sensitivity 
to  speed  30/3U,  the  sensitivity  to  temperature  3Q/3(A9)  and  the  ratio  of 
speed-to-temperature  sensitivity  are  readily  determined.  The  coating  and 
the  substrate  each  reduce  the  quasi-steady  sensitivity  to  speed  and  to 
temperature  oscillations  of  the  ambient  fluid.  The  substrate  and  the 
coating  also  reduce  the  ratio  of  speed-to-temperature  sensitivity  below 
its  maximum  value  of  A9/2U.  The  effect  of  a finite  coating  thickness 
can  be  very  significant  in  water.  At  15  cm/sec  a typical  glass  coating 
of  0.0028  cm  on  a 0.017  cm  diameter  microbead  thermistor  reduces  the 
speed  sensitivity  by  a factor  of  10  and  the  ratio  of  speed-to-temperature 
sensitivity  by  a factor  of  3.2  from  its  value  in  the  absence  of  a coating. 
This  effect  increases  with  increasing  speed.  A typical  glass  coating 
of  I or  2 x IQ-4  cm  over  a thin  film  probe  reduces  the  temperature  of 


- 15  - 


the  wetted  surface  by  approximately  20%.  The  glass  substrate  of  a hot 
film  probe  can  carry  about  20%  of  the  total  heat  flux  from  the  film. 

The  effect  of  the  coating  and  the  substrate  depends  strongly  on  the 
relative  conductivity  of  the  fluid.  There  is  a strong  coating  effect 
but  a moderate  substrate  effect  on  probes  in  water  whereas  in  air  the 
coating  effect  is  not  very  important  but  the  substrate  effect  is 
i mportant. 

A knowledge  of  a probe's  sensitivity  to  both  speed  and 
temperature  is  crucial  for  oceanic  work  because  temperature  activity 
often  accompanies  velocity  activity.  Some  measurements  of  the  intensity 
of  temperature  and  velocity  turbulence  are  now  available.  System  noise 
and  the  probe's  working  life  are  two  other  important  and  sensitivity 
related  considerations.  An  a priori  evaluation  of  a sensor's  perfor- 
mance is  possible.  The  expected  signal  can  be  used  to  judge  a sensor's 
s i gna I -to-noi se  ratio,  its  signal  contamination  and  its  performance 
relative  to  other  sensors.  Because  the  temperature-heat  flux  relation- 
ship is  sensitive  to  geometry  the  merits  of  various  coating  materials, 
coating  thicknesses,  and  probe  shapes  and  sizes  can  be  determined  from 

I 

the  appropriate  function  h(Rg2). 

Oceanic  microscale  temperature  gradients  can,  at  some  times 
and  in  some  locations,  be  small  enough  to  prohibit  their  resolution 
with  the  best  present  technigues  (Gregg,  et  al,  1978  using  thermistors). 
An  increase  of  the  s i gna I -to-noi se  ratio  is  therefore  important. 
Increasing  a thermistor's  current  increases  its  s igna I -to-noi se  ratio 
but  also  increases  its  sensitivity  to  velocity  fluctuations.  Because 


the  heat  flux  relationship 
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= h(R  s) 

Q e 

is  independent  of  the  level  of  self  heatinq  it  applies  to  thermometers 
as  well  as  anemometers.  The  limit  to  which  the  current  through  a 
thermistor  used  for  thermometry  can  be  increased  is  governed  bv  its 
relative  ve I oc i ty-to-temperature  sensitivity,  the  relative  velocity- 
to- temperature  activity  in  the  fluid  being  probed  and  the  level  of 
contamination  of  the  temperature  signal  that  one  is  willing  to  accept. 
Circumstances  may  vary  greatly.  By  using  the  observed  values  for  the 
rates  of  temperature  and  kinetic  energy  dissipation  in  the  ocean  it  is 
shown  in  appendix  C that  temperature  measurements  made  with  glass  coated 
microbead  thermistors  may  be  contaminated  by  velocity  in  the  mixed  layer 
and  in  other  regions  of  high  velocity  but  low  temperature  activity  when 
the  over-heat  A9  is  only  MO-2  °C.  Heated  anemometers,  such  as  grossly 
heated  thin  film  probes  and  thermistors,  are  also  sensitive  to  tempera- 
ture. Several  people  have  suggested  warmed  thermistors  for  high 
frequency  thermometry.  However,  because  even  low  levels  of  internal 
heating  may  cause  velocity  contamination  it  is  not  likely  that  even 
moderately  heated  probes  are  practical  for  thermometry  in  the  ocean. 

As  a final  comment  to  this  section,  calibration  data  is  often 
presented  in  terms  of  a bridge  voltage  and  an  over-heat  ratio.  This 
makes  calibration  comparisons  difficult.  To  determine  the  heat  flux 
versus  flow  relationship  one  then  needs  some  knowledge  of  the  probe's 
electronics  and  the  temperature-res i stance  calibration  of  the  probe. 
Often  the  needed  information  is  not  published.  This  is  acceptable  if 
the  only  purpose  of  publication  is  to  demonstrate  that  a given 


operational  system  works  but  is  not  acceptable  if  the  purpose  is  to 
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demonstrate  the  fundamental  characteristics  of  a probe.  Understand i ng 
of  probe  character i st ics  might  proceed  more  swiftly  if  calibration  data 
were  published  in  terms  of  relevant  parameters  such  as  over-heat  A0 , 
heat  flux,  Reynolds  numbers,  Prandtl  numbers,  and  sensor  dimensions 
instead  of  bridge  voltages  and  over-heat  ratios. 

I 


4.  Unsteady  Heat  Transfer 

The  Nusselt  number  is  a steady  heat  transfer  concept.  The 
response  of  sensors  is  not  governed  by  the  Nusselt  number  when  the 
viscous  boundary  layer  is  unsteady.  An  understanding  of  the  unsteady 
characteristics  of  ideal  probes  can  be  used  advantageously  to  under- 
stand the  unsteady  behaviour  of  real  probes. 


4.1  I dea I Probes 

The  behaviour  of  the  unsteady  viscous  boundary  layer  over  a 
flat  plate  has  been  examined  by  the  author  (appendix  A),  Chenq  and 
Elliott  (1957),  Hill  and  Stenning  (I960),  Lighthill  (1954),  Ostrach 
(1955),  Rott  and  Rosenzweiq  (I960)  and  others.  Only  Ostrach  (1955) 
and  the  author  have  extended  their  examination  to  the  unsteady  heat 
flux  from  a plate.  Ostrach's  results  are  limited  to  very  low  frequencies 
whereas  the  author's  results  apply  to  all  freguencies.  Lighthill  (1954) 
examined  the  unsteady  viscous  boundary  layer  and  the  heat  flux  over 
the  frontal  surface  of  a cylinder  (two-dimensional  stagnation  flow). 

The  author  (appendix  B)  has  examined  the  unsteady  boundary  layer  and 
the  heat  flux  over  a cone  of  half  angle  15°  under  axi-symmetr ic  flow. 

The  boundary  layer  over  a flat  plate  of  length  L departs 
from  its  quasi-steady  behaviour  when  the  frequency  parameter  wL/U>0.3 
where  w is  the  angular  frequency  in  rad i ans/second . Any  criterion  for 
unsteady  behaviour  is  (necessarily)  arbitrary.  The  frequency  parameter 
chosen  is  that  for  which  the  unsteady  surface  shear  stress  on  any  part 
of  the  surface  exceeds  the  quasi -steady  stress  by  more  than  20$.  The 
boundary  layer  over  a cone  of  half  angle  less  than  approximately  20° 
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becomes  unsteady  when  wL/U>0.9  and  becomes  unsteady  over  the  front  of 
a cylinder  of  diameter  D when  u>D/U£5. 

An  examination  of  the  unsteady  first  order  perturbation 
equations  of  an  ideal  constant  temperature  sensor,  like  the  plate 
discussed  in  appendix  A,  is  instructive.  The  equations  are 


3u(  3u | 3uq  3u 

iwui  + uo  i7+  vo  9y“+  ui  aT+  vi  af 


3^  u | 

iwUo  + V 


(10) 


0 at  y = 0 and  x > 0 


u,  = U 


at  y 


3T, 


3T, 


92T . 


i u>T . + u — — + v 

I o 3x  o 3y 


3y 


7 = 


- u. 


3T 

o 

3x 


- v , 


3T 

c 

3y 


(in 


and 


T = 0 at  y = 0,  00  and  x > 0 


. . , 3 9 39  3 ?0  n 

i wA  + u r — + v r—  - < r — y - 0 
o 3x  o 3y  3y^ 


(12) 


9 = 0 at  y = 0,  x > 0 
9 = I at  y * ® 

from  equations  5(a),  16  and  17  in  appendix  A, 

where  (u,v)  = (u  , v ) + eela)+(u,  ,v.  ) 

o o II 

T T , iu>tT  , iu>tn 
T = T + ee  T.  + ye  0 
o i 


U = U (I  + ee'w+) 
o 

is  the  free  stream  speed,  T(x,°°)  = ye'w^  is  the  free  stream  temperature 
and  the  constants  y and  e are  small  compared  to  unity.  The  boundary  layer 
velocity  is  represented  by  the  steady  field  (u^v^)  and  a small  perturba- 
tion ce'^fuj  ,V|  ) due  to  a small  fluctuation  of  the  main  stream  eII^  e'u+. 
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The  steady  temperature  field  due  to  steady  streaming  and  the  constant 

surface  temperature  is  represented  by  T , ee'^T.  is  the  temperature 

field  due  to  the  unsteady  velocity  perturbation,  and  yhe'1'  represents 

the  temperature  field  due  to  a small  independent  perturbation  of  the 

main  stream  temperature  ye'w^.  Because  of  the  non-homogeneous  term 

ioU  in  eguation  10,  which  represents  the  unsteady  pressure  gradient 
o 

in  the  fluid,  the  unsteady  boundary  layer  velocity  ( u ( , v ( ) near  the 
plate  increases  in  magnitude  and  advances  in  phase  from  its  quasi-steady 
value.  A comparison  of  equations  II  and  12  makes  it  quite  clear  that 
the  anemometric  response  (Tj  ) must  be  quite  different  from  the  tempera- 
ture response  (9).  Equation  10  is  non-homogeneous  but  has  homogeneous 
boundary  conditions.  The  heating  source  that  generates  the  perturbation 
temperature  field  T|  is  not  at  the  boundaries  but  is  intrinsically 
distributed  throughout  the  entire  boundary  layer  and  results  from  the 

advection  of  the  steady  boundary  layer  temperature  gradient  (9T  /9x, 

o 

9T  /9y)  by  the  unsteady  boundary  layer  velocity  (u.,v  ).  This  heating 
o i ' 

source  u,*VT  is  a function  of  Prandtl  number,  frequency  and  position. 

- I o 

Proper  frequency  response  calibration  requires  a faithful  simulation 
of  the  forcing  term  (u ,'VT  ) which  can  be  provided  by  vibratinq  the 
probe  longitudinally  in  a steady  flow  (Lighthill,  1954).  Calibrations 
in  a steady  flow  made  by  internally  modulating  the  probe's  temperature 
by,  say,  an  unsteady  electric  current  (Freymuth,  1977,  1978;  Ling,  1955; 
Tan-atichat,  1973;  Gust,  1976  and  others)  can  infer  only  the  quasi- 
steady response  because  the  unsteady  boundary  layer  is  absent  (details 
are  discussed  in  appendix  A). 


MMH 
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The  temperature  response  equation  (II)  is  homogeneous  and 
driven  by  a non-homoqenei ty  at  the  boundary  y = <*>.  This  heating 
source  is  outside  of  the  boundary  layer  and  can  only  penetrate  the 
boundary  layer  by  steady  advection  and  by  diffusion.  Independent  of 
the  level  of  heating,  the  anemometric  response  bandwidth  is  wider 
than  the  temperature  response  bandwidth  because  (I)  the  heating  source 
U|*VT  for  the  anemometric  response  increases  in  magnitude  whereas  the 
heating  source  for  the  temperature  response  is  constant  and  (2)  the 
heating  source  for  the  anemometric  response  is  intrinsic  to  the  boundary 
layer  but  the  heating  source  for  the  temperature  response  is  external 
to  the  boundarv  layer.  The  -3  dB  response  bandwidth  for  a flat  plate 
in  the  Prandtl  number  range  0.7  < Pr  < 15  is 


u)L 

U 


3 dB 


3.7  (I  + 2/Pr) 
(I  + /Pr)2 


for  the  anemometric  response  and 


u>L 

u 


3 dB 


0.56  Pr“-38 


for  the  temperature  response.  Wedges  and  cones  should  have  a wider 
bandwidth  because  their  boundary  layers  are  thinner.  Lighthill  (1954) 
found  that  the  -3  dB  bandwidth  of  the  heat  flux  from  the  front  of  a 
cylinder  in  response  to  speed  oscillations  is  approximately  wD/U  = 5 
at  a Prandtl  of  0.7  (air).  At  a Prandtl  number  of  0.7  the  bandwidth 
of  a cylinder  is  1.7  times  larger  than  the  bandwidth  of  a plate  (if  D=L). 
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The  above  bandwidths  are  independent  of  the  level  of  steady  heatina 
but  apply  strictly  only  to  constant  temperature  surfaces.  For  variable 
temperature  probes  the  thermal  inertia  of  the  probe  may  also  be  impor- 
tant and  further  limit  the  response  bandwidth.  For  either  case,  constant 
temperature  or  variable  temperature  probes,  the  anemometric  bandwidth 
must  be  larger  than  the  temperature  response  bandwidth  because  of  the 
boundary  layer's  behaviour.  However,  the  difference  between  the  two 
responses  will  be  less  pronounced  for  variable  temperature  probes  and 
the  amount  of  this  difference  will  depend  on  the  relative  importance  of 
the  probe's  inertia. 

4.2  Real  probes 

The  high  frequency  heat  flux  from  the  two  conical  probes 
examined  in  appendix  B differs  significantly  from  its  quasi-steady 
value.  The  heated  film  on  these  probes  spans  the  entire  circumference 
at  a position  "downstream"  from  the  apex.  The  film  is  narrow  compared 
to  its  distance  from  the  apex.  These  probes  have  a frequency  response 
which  at  first  seems  surprising;  the  response  rises  with  increasing 
frequency.  Because  of  the  film's  aspect  ratio  the  unsteady  heat  flux 
is  proportional  to  the  viscous  surface  stress  over  the  film  which 
increases  in  magnitude  and  phase  advance  with  increasing  value  of  the 
frequency  parameter  Q = ioL/3U.  Increasing  the  relative  width  of  the 
film  by  either  moving  the  film  towards  the  apex  or  by  widening  the 
film  would  tend  to  diminish  the  rise  of  the  freauency  response.  If 
the  entire  cone  of  15°  half-angle  were  covered  by  the  film  then  the 
frequency  response  should  be  similar  to  the  response  of  a flat  plate. 
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Significant  departure  from  auasi -steady  behaviour  of  the  heat 
flux  from  a wedge  shaped  sensor  in  air  at  frequencies  where  the  boundary 
layer  must  still  be  quasi-steady  have  been  reported  by  Bellhouse  and 
Schultz  (1967).  This  effect  is  caused  by  the  substrate.  The  heat  flux 
into  the  substrate  must,  somewhere,  pass  through  the  probe's  boundary 
layer  to  enter  the  ambient  fluid.  If  the  substrate's  conductivity  is 
large  compared  to  the  fluid's  conductivity,  the  substrate's  heat  flux 
is  sensitive  to  the  flow  because  of  the  boundary  layer  in  its  path.  The 
substrate's  unsteady  heat  flux  in  response  to  speed  oscillations  can 
decrease  with  increasing  frequency  because  of  the  substrate's  finite 
thermal  diffusivity.  If  the  substrate's  conductivity  is  not  large  compared 
to  the  fluid's  conductivity  then  the  heat  flux  through  the  substrate 
is  nearly  independent  of  the  flow  and  does  not  respond  to  speed 
modulations.  The  heat  flux  through  the  substrate  in  response  to 
oscillations  of  the  ambient  fluid's  temperature  on  the  other  hand, 
must  decrease  with  increasing  frequency  regardless  of  the  value  of 
the  substrate- to-f I u i d conductivity  ratio.  Low  frequency  attenuation 
of  the  heat  flux  in  response  to  speed  oscillations  is  evidently  small 
in  water  but  large  in  air.  The  unsteady  calibration  data  of  the  two 
conical  hot  film  probes  discussed  in  appendix  B show  very  little  low 
frequency  attenuation.  Frequency  response  calibrations  in  water  by 
Nowell  (1974)  do  not  show  any  low  frequency  attenuation. 

Because  the  response  bandwidth  is  wider  for  speed  than  for 
temperature  and  because  the  substrate  may  affect  a low  frenuency 
attenuation  of  the  temperature  response  but  not  the  speed  response 
(in  water),  the  ratio  of  temperature-to-veloc ity  sensitivity  is  highest 
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at  zero  frequency.  An  upper  bound  for  the  temperature  contamination 
of  a measured  velocity  spectrum  can  be  calculated  with  the  temperature 
spectrum  and  the  quasi-steady  sensitivity  ratio.  Conversely,  a lower 
bound  for  the  velocity  contamination  of  the  temperature  spectrum  can 
be  calculated  with  the  quasi-steady  value  of  the  sensitivity  ratio. 

The  high  frequency  response  bandwidth  of  thin  film  anemometers 
is  given  approximately  by  oL/U  'v  | and  depends  somewhat  on  the  Prandtl 
number.  If  one  assumes  that  a probe  traverses  through  a "frozen" 
turbulence  field  at  a speed  U then  the  largest  wave  number  o that 
can  be  resolved  because  of  frequency  response  limitations  is 

n ^ - 

0 U 

which  implies  that  oL  ^ I . Thus  the  scale  size  corresponding  to 
frequency  response  limitations  is  also  approximately  the  scale  size 
corresponding  to  spatial  response  limitations. 
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5.  Conclusions 


This  thesis  has  discussed  the  steady  and  unsteady  heat  transfer 
from  various  types  of  sensors.  The  salient  conclusions  are  as  follows. 


5.1  Steady  Heat  Transfer 

(i)  The  Nusselt  number  is  a steady  heat  transfer  concept  and  is 

limited  to  the  frequency  range  where  the  boundary  layer  is  quasi-steady 
A functional  relationship  for  the  steady  heat  transfer  from  probes  with 
an  insulating  coating  and  a substrate  can  be  derived  by  incorporating 
the  concept  of  a Nusselt  number  governed  heat  flux  from  the  wetted 


surface. 


The  functional  relationship  for  the  steady  heat  flux  has 


been  explicitly  derived  for  flake  thermistors,  microbead  thermistors 
and  conical  hot  film  probes.  The  functional  relationship  is  AT/O  = 

ill 

h(R  2)  where  h(R  2)  i (A  + BR  2)_1  because  of  the  thermal  effects  of 
e e e 

the  coating  and  the  substrate. 

I 

(iii)  The  function  h(R  J)  can  be  used  to  calculate  the  quasi- 

e 

steady  sensitivity  of  thermometers  and  anemometers  as  well  as  the 
quasi-steady  contamination  of  thermometers  by  velocity  and  anemometers 
by  temperature. 

(iv)  The  coating  and  the  substrate  each  reduce  the  quasi-steady 
sensitivity  to  speed  and  temperature  oscillations  of  the  ambient  fluid 
as  well  as  the  ratio  of  speed-to- temperature  sensitivity.  The  maximum 
value  of  the  ratio  of  speed-to-temperature  sensitivity  under  forced 
convective  heat  transfer  is  A0/2U. 


(v)  With  some  knowledge  of  h(Rg2)  an  a priori  evaluation  of  a 


i 


\ ' 
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sensor's  performance  is  possible.  The  expected  signal  can  be  used  to 
judge  a sensor's  signa l-to-noi se  ratio,  its  level  of  contamination,  and 

i 

its  merits  relative  to  other  sensors.  Because  the  function  h(R  2)  is 

e 

sensitive  to  aeometry  the  relative  value  of  various  coat i na  materials, 
coating  thicknesses  and  probe  configurations  can  be  assessed.  The 
function  h(Rg?)  is  independent  of  the  probe's  electronics.  Inter- 
comparison of  published  calibration  data  would  be  easier  if  it  were 
presented  in  terms  of  relevant  non-electrical  parameters  such  as  over- 
heat A0,  Reynolds  numbers,  and  Prandtl  numbers  rather  than  bridge 
voltages  and  over-heat  ratios. 

5.2  Unsteady  Heat  Transfer 

(i)  The  viscous  boundary  layer  over  a probe  departs  from  its 
quasi -steady  behaviour  when 

(1)  u)L/U  > 0.3  for  a flat  plate, 

(2)  wL/U  > 0.9  for  cone  of  half-angle  less  than  approximately 
20°,  and 

(3)  u)L/U  > 5 for  a cylinder 

where  L is  the  stream-wise  dimension  of  the  probe  or  the  diameter  of 
the  cy I i nder . 

(ii)  The  heating  source  that  generates  the  unsteady  anemometric  effect 
is,  to  a first  order  perturbation,  the  advection  of  the  steady  boundary 
layer  temperature  gradient  by  the  unsteady  boundary  layer  velocity, 
namely  u,‘VT  . This  heating  source  is  intrinsically  distributed  throuqh- 
out  the  boundary  layer  and  is  a function  of  frequency,  Prandtl  number 
and  position. 


u 
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(iii)  Frequency  response  calibrations  must  realistically  simulate 
the  generation  term  U| -VTq  when  the  boundary  layer  is  unsteady;  probe 
vibration  does  but  internal  heating  does  not. 

(iv)  The  response  of  the  heat  flux  to  oscillation  of  the  ambient 
fluid's  temperature  is  generated  by  the  unsteady  temperature  boundary 
condition  at  the  exterior  of  the  boundary  layer.  Consequently  the 
response  bandwidth  is  larger  for  speed  oscillations  than  for  temperature 
oscillations.  Both  bandwidths  depend  strongly  on  the  Prandtl  number. 
Because  both  bandwidths  are  independent  of  the  level  of  internal  heating 
the  above  conclusion  also  holds  for  nearly  unheated  temperature  sensors. 

(v)  The  heat  flux  into  the  substrate  is  sensitive  to  the  boundary 
layer  if  the  ratio  of  substrate- to-f I u i d conductivity  k /k^  is  large. 

If  k /k,  is  not  larqe  then  the  heat  flux  into  the  substrate  is  nearly 
independent  of  the  flow.  If  k^/k^  >>  I then  the  unsteady  heat  flux 
into  the  substrate  in  response  to  speed  oscillations  decreases  with 
increasing  frequency.  Regardless  of  the  value  of  k^/k^.  the  heat  flux 
into  the  substrate  in  response  to  temperature  oscillations  of  the 
ambient  fluid  decreases  with  increasing  frequency. 

(vi)  The  ratio  of  temperature- to-ve loc i ty  sensitivity  is  highest 
at  zero  frequency. 

(vii)  An  upper  bound  for  the  temperature  contamination  of  a measured 
velocity  spectrum  can  be  calculated  with  the  temperature  spectrum  and  the 
quasi-steady  sensitivity  ratio.  Conversely,  a lower  bound  for  the 
velocity  contamination  of  the  temperature  spectrum  can  be  calculated 
with  the  quasi-steady  value  of  the  sensitivity  ratio. 
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Abstract 

The  unsteady  viscous  stress  on  a flat  plate  computed  numerically 
from  a simplified  unsteady  momentum  equation  that  is  not  restricted  by 
frequency  compares  favourably  with  the  frequency  restricted  calculations 
available.  The  unsteady  heat  flux  from  the  same  plate  maintained  at  a 
constant  temperature  was  also  computed.  The  heat  flux  responds  differently 
to  oscillations  of  the  ambient  fluid's  speed  than  it  does  to  oscillations 
of  the  ambient  fluid's  temperature.  The  response  to  speed  oscillations 
has  a maximum  value  at  low  frequencies  and  a larger  bandwidth  than  the 
response  to  temperature  oscillations.  An  examination  of  the  unsteady 
temperature  equations  shows  that  an  oscillation  of  the  plate's  surface 
temperature  does  not  generate  a heat  flux  that  is  comparable  to  the  heat 
flux  in  response  to  either  speed  oscillations  or  to  temperature  oscillations 
of  the  ambient  fluid  when  the  boundary  layer  is  unsteady. 
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I . I ntroduct ion 

A significant  motive  for  understanding  the  unsteady  heat  flux 
from  bodies  exposed  to  streaming  comes  from  the  common  use  of  heated  sensors 
for  high  frequency  turbulence  flow  measurements.  Spectra  of  velocity  and 
temperature  may  be  significantly  biased  (filtered)  if  the  sensor's  un- 
steady response  is  different  from  its  quasi-steady  response,  in  which 
case,  some  data  correction  is  required.  Wire-shaped  sensors  have  tradition- 
ally been  used  in  air.  In  water,  however,  it  is  more  common  to  use  thin 
film  sensors  that  are  wedge-shaped  (Nowell,  1974),  cone-shaped  (Gargett, 
1976,  Grant  et  al.,  1962,  1968),  and  cylindrical. 

The  dimensions  of  film  type  sensors  differ  by  several  factors  of 
10  from  wire  sensors;  wires  have  typical  diameters  of  ICT6  meters,  cylin- 
drical film  sensors  I 0-4  meters,  and  wedge  and  cone  shaped  sensors  have 
films  about  JO-3  meters  long.  As  shown  later,  a frequency  parameter  that 
characterizes  the  unsteady  boundary  layer  is  coL/U^  where  w is  the  angular 
frequency  of  oscillation  (in  radians  per  second),  L is  a relevant  dimension 

(diameter  or  lenqth  in  the  flow  direction),  and  U is  the  ambient  fluid's 

o 

steady  speed.  If  wL/U^  <<  I the  flow  is  essentially  quasi-steady,  its 
behaviour  at  any  instant  is  that  appropriate  to  steady  conditions  at  the 
instantaneous  value  of  the  fluid's  speed  and  temperature.  If  wL/Uo  2 I 
the  flow  is  unsteady  and  may  behave  differently  from  the  quasi-steady  flow. 
It  is  readily  seen  that  wire  type  sensors  operate  in  the  quasi-steady 
frequency  range  for  all  practical  purposes.  However,  the  boundary  layer 
on  cylindrical  films  is  unsteady  when  w z I 04  Uq  (Uo  in  meters/second ) and 
is  unsteady  on  wedge  and  cone  shaped  sensors  when  w 2 I 03  Uq.  Thus  unsteady 
effects  are  within  the  frequency  domain  of  typical  measurements  made  with 
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wedge  and  cone  shaped  film  turbulence  sensors.  Unsteady  boundary  layer 

effects  have  been  reported  by  Nowell  (1974)  and  Evans  (1963).  In  both 

cases  the  response  to  external  speed  oscillations  was  measured  by  vibrating 

the  probe  in  the  direction  of  the  steady  stream.  Nowell  showed  that  the 

unsteady  response  of  wedge  shaped  sensors  diminished  with  increasing 

frequency  when  w/U  s IO-3metre  1 . The  response  was  down  by  3 dB  at 
o 

120  cyles/sec  for  U =0.31  m/sec  and  500  cycles/sec  for  U =1.12  m/sec. 

o o 

For  cone  shaped  sensors  with  a narrow  film  at  a distance  L from  the  rear 
of  the  apex,  Evans  reported  that  the  response  rises  when  wL/U^  ;>  0.9 

The  frequency  response  of  heated  sensors  is  often  inferred  from 
internal  heating  methods;  the  probe  is  exposed  to  a steady  streaming,  its 
temperature  is  oscillated  by  an  electric  current  and  the  resulting  heat 
flux  is  measured.  Whereas  one  wants  to  know  the  response  to  oscillations 
from  outside  of  the  boundary  layer  one  is  calibrating  the  response  to  an 
internal  oscillation.  The  calibrations  made  by  internal  heating  are  related 
to  the  response  to  external  speed  and  temperature  oscillations  by  implicitly 
assuming  that  the  response  (to  external  oscillations)  is  quasi-steady 
(Freymuth,  1978).  When  the  boundary  layer  is  quasi-steady  it  does  not 
matter  whether  the  oscillations  are  external  or  internal  because  the  inertia 
and  the  thermal  inertia  of  the  boundary  layer  are  negligible.  However,  when 
the  boundary  layer  is  unsteady  the  heat  flux  in  response  to  external 
oscillations  is  different  from  the  heat  flux  in  response  to  internal 
temperature  oscillations;  at  unsteady  frequencies  (uL/U  j I ) the  former 
decreases  with  increasing  frequency  while  the  latter  increases  with  increas- 
ing frequency  (as  shown  later).  At  high  frequencies  the  flux  in  response  to 
internal  oscillations  depends  only  on  the  fluid's  thermal  properties  and  ui, 
whereas  the  flux  in  response  to  external  oscillations  always  depends  on  the 
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fluid's  speed  and  Prand+I  number. 

The  point  here  is  that  if  the  response  of  a sensor  is  not  affected 
by  the  inertia  or  the  thermal  inertia  of  the  boundary  layer  then  internal 
heating  can  give  useful  information.  If  the  sensor  is  operated  in  a 
constant  temperature  mode  then  internal  heating  can  tell  about  the 
electronics.  If  the  sensor  is  operated  in  a constant  current  mode  then 
internal  heating  can  measure  the  thermal  inertia  of  the  probe  and  the 
limitation  of  the  electronics.  But  if  the  response  is  limited  by 
boundary  layer  processes  then  internal  heating  cannot  give  any  infor- 
mation about  the  effect  of  the  boundary  layer  because  (I)  the  unsteady 
velocity  in  the  boundary  layer  is  not  present  during  calibration  and 
because  (2)  the  heat  flux  from  the  probe  responds  differently  to  external 
oscillation  than  to  internal  oscillation.  The  probe's  constant  temperature 
electronics  cannot  control  the  behaviour  of  the  boundary  layer.  Hence, 
frequency  response  calibrations  of  cylinder,  wedge  and  cone  shaped  sensors 
/ made  by  internal  heating  may  be  erroneous  in  the  frequency  range  where 

the  boundary  layer  is  unsteady. 

In  practice  the  response  may  be  complicated  by  substrate  effects 
(Bellhouse  and  Schultz,  1967),  coating  effects  and  possibly  limitations  of 
the  electronics.  However,  an  examination  of  the  inherent  hydrodynamic 
limitations  of  an  ideal  sensor  like  a flat  plate  has  practical  implications 
because  these  limitations  must  also  apply  to  real  sensors  and  other  heated 
surfaces  exposed  to  an  unsteady  streaming. 

A considerable  amount  of  work  has  been  published  on  the  plate's 
non-steady  velocity  boundary  layer  (Cheng  and  Elliott  1957,  Hill  and  Stenning 
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I960,  Rott  and  Resenzweig  I960).  However,  little  of  this  work  has  been 
extended  to  the  unsteady  response  of  the  heat  flux  to  speed  or  to  temperature 
oscillations  in  the  ambient  stream.  Liqhthill  (1954)  investigated  the 
unsteady  flow  over  a flat  plate  for  the  case  of  small  oscillations  in  the 
magnitude,  but  not  the  direction,  of  the  ambient  fluid.  By  using  the 
Mrm^n-Poh I hausen  method,  Liqhthill  derived  an  approximate  analytic  solution 

I 

l t 

for  the  velocity  field  for  low  frequencies.  The  unsteady  heat  flux,  however, 

could  not  be  derived  by  the  same  method,  mainly  because  of  the  plate's  zero 

steady  pressure  gradient.  Ostrach  (1955),  using  a numerical  method  valid 

for  low  frequencies  (wL/U  <<  I),  found  that  the  heat  flux  increases  slightly 

o 

in  magnitude  and  has  a phase  advance  over  its  quasi-steady  value. 

The  author  is  unaware  of  any  work  on  the  response  of  heated  sensors 
to  temperature  oscillations  of  the  ambient  fluid. 

In  the  following,  attention  is  restricted  to  the  case  where  the 
magnitude  but  not  the  direction  of  the  ambient  fluid  oscillates  about  a 
steady  mean,  and  where  the  amplitude  of  oscillation  is  small  enough  to 
avoid  harmonic  distortions.  The  Reynolds  number  range  is  restricted  to 
where  the  boundary  layer  approximations  are  valid  and  where  the  boundary 
layer  is  laminar.  It  is  assumed  that  the  temperature  field  does  not 
influence  the  velocity  field.  The  surface  temperature  of  the  plate  is 
taken  as  constant  at  0 and  the  amplitude  of  oscillation  of  the  ambient 

fluid's  temperature  is  assumed  to  be  small  compared  to  0.  The  analysis 

] l 

has  applications  for  constant  temperature  sensors  and  other  surfaces  whose 
temperature  remains  constant  or  practically  constant  when  exposed  to 
unsteady  streaming. 


The  viscous  boundary  layer  is  examined  in  section  2.  Section  2a 


deals  with  the  steady  boundary  layer  velocity  and  2b  deals  with  unsteady 
flow  in  the  boundary  layer.  The  heat  flux  is  calculated  in  section  3. 

The  unsteady  flux  in  response  to  speed  oscillations  in  ambient  fluid  is 
calculated  in  3a  and  the  response  to  ambient  fluid  temperature  oscillations 
is  examined  in  3b. 


! 
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2.  The  Viscous  Boundary  Layer 


The  boundary  layer  equations  for  two  dimensional  incompressible 


flow  with  an  ambient  fluid  speed  U(t)  paralle 

1 to  the  plate  are  (Lighthill 

1954) 

9u  u9u  + v3u  _[_  9£  _ 92u 

9t  9x  9y  p 3x  9 y^ 

(la) 

^ + & = 0 

9x  9y 

(lb) 

— + 1 — = 0 

9t  p 9x 

(Ic) 

with  boundary  conditions 

u=v=0@y=0  for  x > 0 and 

u •>  U(t)  as  y + “> 

where  p is  the  fluid  density  and  v the  kinematic  viscosity  (Figure  1). 

We  now  let  the  quantities  U,  u,  and  v perform 

small  oscillations  about  a 

steady  mean  and  write 

U(t)  = U (l+eelll)+) 

o 

(2a) 

u = uQ(x,y)  + eUj  (x,y,u))eiu>+ 

(2b) 

v = vQ(x,y)  + ev | (x,y,w)e'w^ 

(2c) 

where  e <<  I is  the  relative  amplitude  of  the  speed  oscillations  and  is 


therefore  a suitable  perturbation  parameter. 

2a . Steady  Flow 

The  well  known  Blasius  transformation  given  in  Schlichting  (1955, 
p.  103),  reduces  the  steady  momentum  equation  (la)  to 

I 1 Ml 

ff  + 2f  = 0 
f (0)  = t' (0)  = 0 

f ' (»)  = I , (3) 


4 
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where 


vx 


u = U f (n ) , 
o o 


and 


U 


v0(x,y) 


vx 


U 


(nf 


f ) 


(4) 


and  primes  denote  differentiation  with  respect  to  n . 

2b.  Non-Steady  Flow  e « I 

By  substituting  equations  (2a,b,c)  into  (la,b,c)  and  by  only 


retaining  terms  up  to  order  e one  gets  (Lighthill,  1954) 


3 u , 


3u , 


3u 


3 u 


iuu.  + u — — L + v — - + u.  t — ' + v.  r— ^ = iuU  + v 
I o 3x  o3y  I 3x  I 3y  o 


32u, 


3y^ 


(5a) 


3 u . 


3v, 


3x 


3y 


= o 


(5b) 


where  u.  = v = 0 @ y = 0 and  x > 0 and 


u . = U as  y -*■ 
o 


Lighthill  (1954)  found  an  approximate  solution  to  equation  (5a)  for  small 
values  of  the  frequency  parameter  9 = tox/U  . For  very  large  values  of  the 


parameter  tox/U^  Lighthill  reduced  equation  (5a)  to  a "shear  wave"  equation 


32  u . 


iw(u . - U ) = v — — 5- 
I o 3yl 


(6) 


which  has  the  solution 


u,  = U { I 
o 


, I iu.-,  wx 

exp(-yv— )} , — » 

o 


(7) 


The  substitution  of 

1 

u. (x,y,u)  = U g (n,0) 


(8a) 


U 


V|  (xfy,u)  = ^ /jp  (ng  - g - 29 
0 


la} 


3n 


(8b) 
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which  satisfies  continuity  (5b),  into  (5a)  yields 


where 


1 1 i | 1 1 t | 1 1 

g + j fg  - ing  + 2 f g - 0 = -in 

% » f I f 

0 = n (f  9 - f a5 


g(0,R)  = g (0,R)  = 0 

g((“»n)  = I 

g(n,0)  = -L  (f  + nf  ) • 


Rott  and  Rosensweig  (I960)  have  partially  solved  equation  (9b)  numerically 
by  writing  g as  a power  series  in  SI,  that  is 


g(n,R)  = £ (iR)nh  (n) 

n=0  n 


for  R < I , 


and  calculated  h^  and  h|.  However,  a power  series  expansion  of  the  above 
type  is  impractical  for  R n,  | . Cheng  and  Elliott  (1957)  have  solved  the 

momentum  equation  (equation  (5a))  to  higher  orders  in  R and  have  tabulated 

1 1 

sufficient  coefficients  to  evaluate  the  dimensionless  surface  stress  g (0,0) 
to  the  fourth  order  in  R . 

Equation  (9a)  becomes  an  ordinary  differential  equation  where  R 
is  only  a parameter  if  the  term  0 in  equation  (9a)  can  be  ignored.  Clearly 
as  R ->  0 0 -*■  0.  For  very  high  frequencies  (R>>l)  we  have  from  equation  (7), 
(8a),  and  (9b)  that 

i t 1 1 1 1 t 1 1 

0 = — (nf  g - nf  q + f g).  (10) 

For  n < 2 the  function  f i s to  a good  approximation  (Schl ichting,  1954,  p.  107) 
i 1 1 

f = j n2f  (0) 

which  makes  the  term  0 for  n<2 

0 ^ 9 f (0)  (n2g  - ng  + g) 


I 
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whereas  the  advective  term  is 


t » it 


{ fg  + f g}  % 7 f i°)  n2g  + g} . 


Because  the  magnitude  of  g is  never  much  larger  than  unity  the  term  0 
is  comparable  to  the  advective  terms  for  n<2  at  hiah  frequencies.  For 
n<4  an  adequate  approximation  for  f is  (Schlichting  1954,  p.  107) 
f = n 

which  makes  the  term  0 approximately 

0 £ j ng' ' 

and  equal  to  the  advective  term.  However,  for  high  frequencies  the  ad- 
vective term  is  small  relative  to  the  pressure  gradient  in,  the  inertia 

i mi 

term  -iQg  and  the  diffusion  term  a . The  resulting  equation,  valid 


for  fi» I , is 


= inCg'-l ) 


which  is  equivalent  to  the  "shear  wave"  equation  6. 

For  low  frequencies,  that  is  ficcl,  a first  order  expression  for 
g is  given  by  Lighthill  (1954)  and  can  be  written  as 

g(n,n)  = g(n,0)  + ifthtn),  where  h is  a real  function  of  order 


unity.  Consequently. 


0 = if)  { f ’ h’  - f”h} 


(I  la) 


whereas  the  advective  term  is 


■=■  { [fg’  ’ (n ,0)  + f”g(n,0)]  + ilMfh"  + f"h)}. 


(I  lb) 


i 


The  term  0 is  small  compared  to  the  advective  term  because  fl<<l  and 
because  the  negative  sign  in  (I  la)  causes  a partial  cancelling  except  for 
large  n where  h’  may  be  slightly  negative.  In  summary  then,  0 is  identically 
zero  at  n=0,  is  small  compared  to  the  advective  terms  for  small  n,  and  is 
comparable  to  the  advective  terms  for  large  Q but  small  compared  to  the 
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dominant  terms  which  are  diffusion,  inertia  and  pressure. 

If  the  velocity  field  does  not  have  any  singularities  with 
respect  to  R,  such  as  resonances  and  notches,  then  one  can  reasonably 
expect  that  the  magnitude  of  0 transits  smoothly,  if  not  monoton ica I I y, 
from  zero  at  R=0  to  a value  comparable  to  the  advective  term  at  very  large 
values  of  R.  Thus  the  deletion  of  the  term  0 from  equation  (9)  cannot 
lead  to  serious  errors  in  computing  g,  its  inclusion  on  the  other  hand 
would  lead  to  formidable  computational  difficulties  near  R=l.  (In 
comparison  the  advective  terms  are  easily  handled  by  numerical  routines.) 

The  simplified  unsteady  momentum  equation 

g* ' ' + Jj-  fg1'  - i Qg  ’ + f''g  = iR 

g(0,R ) = g'(0,R)  = 0,  g'(“,R)  = I (12) 

and  the  steady  momentum  equation  (3)  have  been  solved  numerically  to  an 
accuracy  of  I part  in  I 06  or  better,  by  a collocation  routine  currently 
being  phased  into  operation  at  the  Computing  Science  Department  of  The 
University  of  British  Columbia  (Ascher,  Christiansen,  and  Russell,  1977). 
The  frequencies  R examined  range  from  0 to  82  in  steps  of  (O.l)2.  The 
computed  values  for  f,  f',  and  fM  agree  with  the  values  tabulated  in 
Schlichting  to  within  one  unit  in  the  least  significant  digit.  An  estimate 
of  0 at  R=  I based  on  the  solution  to  equation  (12)  will  be  given  later. 

The  unsteady  velocity  profile  g'  is  plotted  in  Figure  2 for  R=0,  0.25, 

I . 00  and  5.29. 

Figure  3 compares  L i g ht h ill’s  low  frequency  solution  at  his 

critical  frequency  of  R=0.588,  which  in  the  notation  used  here  is 

I - (I  - 0. 1 7ln )2{ ( I - 0. 1 7 1 n - 0.0882n2)  -i  0.5l5n(l  - 0.229n)} 

against  the  computed  velocity  profile  at  the  same  frequency.  The  frequency 


Figure  2.  The  real  and  imaginary  components  of  the  boundary  layer 

perturbation  velocity  profiles  at  values  of  the  frequency 
parameters  fl=0.0,  0.25,  1.0  and  5.29.  The  imaginary 
component  of  g'  is  zero  for  Q=0. 


Figure  3.  A comparison  between  the  computed  perturbation  velocity 
profile  and  L i gh t h ill’s  profile  at  the  "critical" 
frequency  Q=0.588. 
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0=0.588  is  where  Lighthill  joined  his  low  frequency  solution  to  the  shear- 
wave  solution  (7).  The  real  part  of  g'  compares  quite  favourably  with 
L i ghth ill’s  solution,  the  discrepancy  between  the  two  curves  being  approxi- 
mately 10$.  The  imaginary  component  of  g' , however,  differs  significantly 
from  Lighthill's.  Using  the  computed  data  for  g at  <1=1.0  and  0=1.24  it  is 
possible  to  estimate  the  neglected  term  0 and  compare  it  to  the  other  terms 
in  the  momentum  equation,  table  2.  The  magnitude  of  0 is  less  than  1$  of  the 
pressure  gradient  iO  at  n=0.5  and  n=5  and  approximately  5$  at  n=2. 

The  local  steady  stress  on  a plate  is  (Sch I icht i ng) 

AT 

t (x)  = pU  /—  f " (0). 

0 o vx 

Similarly,  the  local  unsteady  stress  is 

nr 

t, (x,n)  = yU  ey  — g' ' (0,0) 

1 o v vx  a 

The  stress  ratio  g' ' (0,0 )/f ' ' (0)e  is  plotted  in  Figure  4(a),  and  included 
for  comparison  is  Lighthill's  low  frequency  formula 
| (I  + 1.7  iO) 

and  the  fourth  order  formula  of  Cheng  and  Elliott 

| {(I  + 0.9429402  - 0.54I0504)  + i ( I . 703530  - 0.7382803)} 

The  magnitude  of  Lighthill's  first  order  stress  formula  is  smaller  than 
the  computed  stress  for  all  frequencies  below  0=0.9  and  the  phase  is 
consistently  larger.  On  the  other  hand  the  agreement  between  Cheng  and 
Elliott's  fourth  order  formula  and  the  computed  stress  is  remarkable, 
the  maximum  difference  amounting  to  less  than  1.6$  in  magnitude  and  1.4$ 
in  phase  for  0<0.64.  For  0>0.7  the  agreement  is  quite  poor;  as  expected 
their  series  expression  diverges.  Table  I summarizes  the  unsteady 


Author 

Chen  C Elliot 

L 1 ghth i 1 1 

fi 

Magnitude 

Phase 

Magnitude 

Phase 

Magni tude 

Phase 

0.0 

1.500 

0.00 

1.500 

0.00 

1.500 

0.00 

0.01 

1.500 

1.02 

1.500 

0.98 

1.500 

0.97 

0.0<t 

1.508 

4.06 

1.506 

3.98 

1.503 

3.89 

0.09 

1.537 

8.90 

1.529 

8.62 

1.517 

8.70 

0.16 

1.607 

14.85 

1.588 

14.75 

1.554 

15.22 

0.25 

1.730 

20.89 

1.703 

21.41 

1.630 

23.00 

0.36 

1.900 

26.24 

1.882 

27.47 

1.759 

31.47 

0.49 

2. 107 

30.69 

2.115 

32.04 

1.952 

39.79 

0.64 

2.345 

34.27 

2.363 

34.69 

2.217 

(47.41) 

0.81 

2.609 

37-06 

2.552 

35.47 

1.00 

2.893 

39-16 

(2.553) 

(34.55) 

1.21 

3-150 

40.69 

1.44 

3.496 

41.78 

1.96 

4.116 

43.11 

3.24 

5.354 

44.19 

9-00 

9.010 

44.84 

64.00 

24.09 

44.99 

Table  Al  The  unsteady-to-steady  stress  ratio  compared  to 

Cheng  6 Elliot's  (1957)  fourth  order  formula  and  L I ghth Ill's 
(195*0  first  order  formula. 
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dimensionless  surface  stress  ratios.  For  very  large  frequencies  the 
non-dimensional  surface  stress  approaches 
g' ' (o,a)  = /7n 

in  agreement  with  the  "shear  wave"  solution  given  in  equation  (7). 

The  total  steady  drag  per  unit  width  on  a plate  of  length  L 
wetted  on  both  sides  is  (Sch I icht i ng,  1955,  p.  108) 

/uT 

D = 4 uU  / — f " (0) . 

0 O V 

Similarly  the  total  unsteady  drag  is 


/U 

D.  = 2eyU  / 

1 0 0 vx 


g' ' (0,0)dx 


/U L | C 

= 4epU  / -2-  ■—  f g’'(0,n)dy^ 


V C*  0 


where  £2  = wL/U  is  the  non-dimensional  plate  frequency.  The  ratio  of 


these  two  drags,  namely 


Us 

eD 


5 


f q'  ' (0,n)  d/7 


o f'MOir  0 

is  plotted  in  figure  4(b),  and  exhibits  characteristics  similar  to  the 
local  stress  ratio,  namely  an  asymptotic  trend  to  /ilf  at  large  frequencies. 


The  uns+eady-to-steady  +ota 
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3.  The  Thermal  Boundary  Layer 

The  heat  flux  from  a plate,  maintained  at  a constant  temperature 
0,  can  be  modulated  by  a change  of  the  ambient  fluid's  speed  and  by  a change 
of  the  ambient  fluid's  temperature.  For  both  cases  it  is  convenient  to 
measure  temperature  with  reference  to  the  ambient  fluid's  steady  tempera- 
ture. The  equations  for  the  distribution  of  the  temperature  T in  a boundary 
layer  without  frictional  heating  are  (Lighthill,  1954) 


31 

3t 


+ u 


II 

3x 


+ v 


3T 

3y 


= K 


ail 

2 

3y 


(13) 


T = 0@y  = Ox>  0 

-r  -r  / t \ i Cot 

I = T (t)  - ye  as  y + ® 

where  K is  the  thermal  diffusivity  of  the  fluid  which  is  taken  to  be  uniform. 
Because  equation  (13)  is  I inear  in  T,  we  can  let  T be  the  sum  of  three 
component  terms.  That  is,  let 

T = TQ(x,y)  + eT|  (x,y,u>)  e'w^  + Y0(x,y,w)  e'u^  (14) 

where  T^  is  the  temperature  distribution  due  to  steady  streaming  over  the 
constant  temperature  plate,  eT  ^ is  the  perturbation  temperature  field  due 
to  the  speed  modulation  ell^  at  frequency  w,  and  y9  is  the  temperature  field 
due  to  an  independent  infinitesimal  modulation  of  amplitude  y of  the  ambient 
temperature  at  the  same  frequency. 

The  substitution  of  equation  (14)  and  the  velocity  expressions 
(2b)  and  (2c)  into  equation  (13)  yields,  to  first  order  in  e,  the  following 
three  sets  of  equations. 


* The  reader  should  be  careful  to  distinguish  0 from  0 in  this  section. 
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Tj  = 0 @ y = 0 
T.  = 0 as  y •+•  « 


iu)0  + u + v - K = 0 (17) 

° 3x  o 3y  ay2 

6 = 0 @ y = 0 

0 = I as  y + “ 

A comparison  of  equations  (16)  and  (17)  makes  it  quite  clear 

that  the  anemometric  response  (T|)  must  be  very  different  from  the 

temperature  response  (0).  Equation  (16)  is  non-homogeneous  but  has 

homogeneous  boundary  conditions.  The  unsteady  heating  source  that  generates 

the  temperature  field  Tj  is  not  at  the  boundaries  but  is  intrinsically 

distributed  throughout  the  entire  boundary  layer  and  results  from  the 

advection  of  the  steady  boundary  layer  temperature  gradient  OT^/Sx, 

3Tq / 3 y ) by  the  unsteady  boundary  layer  velocity  ( u j , v j ) . This  heating 

source  U|.VTq,  as  shown  later,  is  a function  of  Prandtl  number,  frequency 

and  position.  The  temperature  response  (0)  equation  (17),  on  the  other 

hand,  is  a homogeneous  equation  driven  by  a non- homogene i ty  at  the  boundary 

y=“>.  The  heating  source  is  outside  of  the  boundary  layer  and  can  only 

penetrate  the  boundary  layer  through  advection  by  the  steady  velocity 

field  (u  and  v ) and  by  diffusion, 
o o 
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Proper  frequency  response  calibration  requires  a faithful 
simulation  of  the  generating  term,  U|-7T  when  determining  the  anemo- 
metric  response  and  a simulation  of  the  boundary  condition  9=1  at  y^*=° 
when  determining  the  temperature  response.  The  equations  for  an  internal 
temperature  oscillation  of  exp(iwt)  are 

3T.  3T.  32T. 

iuT.  + u - — - + v - — - - K — = 0 
i o 3x  o 3y  3yz 

T / \ i (18) 

T.(x,0,uj)  = I 

T . (x ) = 0. 

For  quasi-steady  frequencies  (wx/Uo<<l)  the  thermal  inertia  of  fhe  boundary 
layer  iwT.  is  small  compared  to  the  other  terms  and  it  is  clear  that 
T.  = 1-0  and,  using  equation  (15),  that 


Using  equation  (16)  one  has  for  ox/Uo<< 


and  hence  there  is  a simple  relationship  between  T ^ , 0,  and  T.  when 
wx/U^cd.  When  aix/U^j  I there  does  not  exist  any  simple  relationship 
between  Tj  , 0 and  T.,  because  (I)  the  generating  term  Uj’VT  becomes 
unsteady  but  is  not  present  in  the  equation  for  T.  and  (2)  the  thermal 
inertia  of  the  boundary  layer  becomes  important.  For  very  high  frequencies 
(ox/Uo>>l)  equation  (18)  reduces  to 

32T . Tin  , . 

.-r  i ^ T(0,O))-l 

i wT . - K — -x—  =0,  T/  , _ 

i 3yz  T(®,w)  = 0 


because  only  terms  involving  w and  the  derivatives  of  highest  order  need 
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+o  be  retained.  The  solution  is 

Tj (y,w)  = exp  {-y  J — } 

The  effect  of  the  internal  temperature  oscillations  on  the  fluid  is 

restricted  to  a layer  whose  thickness  decreases  with  increasing 

1 

frequency.  The  surface  heat  flux  is  k(iw/K)2  which  increases  with 

increasing  frequency  and  does  not  depend  on  the  fluid's  speed  or  Prandtl 

number.  Solutions  for  Tj  and  8 (section  3 a & b)  show  that  the  heat  flux 

in  response  to  speed  or  temperature  oscillations  in  the  main  stream 

decreases  with  increasing  frequency  ..  en  wx/Uo>0.3  and  always  depends 

on  the  speed  and  Prandtl  number  of  the  main  stream.  In  summary,  the 

response  as  measured  by  an  internal  heating  method  can  only  be  related 

to  the  response  to  external  oscillations  when  ox/U  <<l.  When  ux/U  2 I 

0 o 

internal  heating  methods  do  not  correctly  measure  the  response  to  either 
speed  or  temperature  oscillations  in  the  main  stream. 

Using  the  Blasius  transformation,  equation  (15)  can  be  written 
as  ( Sch I icht i ng,  1954,  p.  264) 

T " + Jr  PrfT  ' = 0 
0^0 

T = 0 n = 0 
o 

r - n (l9 

T = 0 n -*■  00 

0 

which  has  the  solution 

00 

To(n,Pr)=-^y  / [f»’(z)]Prdz 
z=n 


where 


Pr  = v/K 
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F(Pr)  / [f  " ( z ) J ' d z 
0 

By  making  use  of  the  B I i i u . trar>:  rmati  >n  and  equation'  (8a),  (8b)  and 
(20),  the  per+ur  ati  n temperature  equation  (K 


- — J ' ' + — fT  ' 
Pr  I 2 1 


T i ( 0,0 ) = 0 

t,  (®,n)  - o 


a lt"<nj'T  (0 , „ ja, , 

2 F 9 ^ 30 


Vn’0)  = U0  = " 2 F [f"(n)] 
0 


and  similarly  the  temperature  response  equation  becomes 


P^9"  + - te'  - iQQ  = 


e(0, n)  = o 
e(°,a)  = i 


e(n,0)  = I 


J f tf”(z)]  dz 
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3a . The  Anemometric  Response 

A perturbation  temperature  field  solution  (T.)  for  very  high 
frequencies  that  is  valid  where  the  steady  temperature  T , is  a linear 
function  of  n was  given  by  L i ght h i i i (1954)  and  is 

Tl  (n,Q)  = in  (nd-PD2  + PrtnN-Pr)  + Me"n'' 

/ift 


2 Pr  -n/iQPr  ^ 
e ) 
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Thus  for  high  frequencies  the  last  term  in  equation  20 

a - 1 >'  <5T,  - "V> 

and  has  the  same  magnitude  as  the  advective  term  ^ f T j ' but  is  small 
compared  to  the  dominant  terms  diffusion,  thermal  inertia  and  heating. 

The  reasons  given  for  dropping  the  term  0 from  the  momentum  equation 
apply  equally  to  the  term  f'Q3T|/3ft  in  equation  (20). 

The  simplified  perturbation  temperature  equation  is 

V +5,PrV  - laPrTi  = 2f^nSTllJ-  '9  + 21  ft > (22> 

Tj  = 0 @ n = 0,». 

The  term  g + 2ft  3g/3ft  can  be  evaluated  from  the  solution  of  the  non-steady 
momentum  equation.  Equation  (22)  has  been  solved  numerically  by  a finite 
difference  formulation  with  881  points  for  n in  the  range  0 to  8.8.  The 
tri-diagonal  matrix  representing  the  finite  difference  term  of  (22)  was 
/ solved  by  Gaussian  elimination.  The  generating  term 

{9+20^1 

which  represents  the  heat  source  for  the  anemometric  effect  is  calculated 
from  the  solutions  to  the  momentum  equations  (4  and  12)  and  is  plotted  in 
Figure  5 for  Prandtl  numbers  I and  8 at  frequencies  0=  0 and  I.  The 
forcing  increases  in  magnitude  and  phase  advance  with  increasing  frequency; 
it  has  vanishing  boundary  values  and  an  amplitude  maximum  with  a value  and 
location  that  depends  strongly  on  the  Prandtl  number.  When  ux/U^jl  the 
unsteady  anemometer  response  can  only  be  measured  by  a method  that 


simulates  this  generating  term. 

The  steady  local  heat  flux  from  the  plate  is 

AT 

-k  / — T ' (0,Pr) 

VX  0 

where  k is  the  fluid's  conductivity.  Similarly  the  local  perturbation 
heat  flux,  that  produces  the  non-steady  anemometric  effect,  is 


/u" 

-ek  / — T.  ' (n,Pr) 


vx 


n=0 


The  perturbation-to-steady  local  heat  flux  ratio 

T ' (f2,Pr)  | 

r=0 


eT  ' (Pr) 
o 


n=0 


which  is  effectively  the  local  unsteady  Nusselt  number  is  plotted  in 
Figure  6(a).  The  heat  flux  magnitude  rises  above  the  quasi-steady  value 
of  0.5  for  all  Prandtl  numbers  from  0.7  to  15.  The  highest  heat  flux 
peak  is  attained  at  a Prandtl  number  of  two,  for  larger  Prandtl  numbers 
the  peak  height  decreases  and  its  location  shifts  to  lower  frequencies. 

A more  detailed  plot  is  in  Figure  6(b).  A power  series  expression  for 
the  local  perturbation  heat  flux  would  have  to  be  of  third  or  higher 
order  to  duplicate  the  low  equency  curves  shown  in  Figure  6(b).  For 
high  frequencies  the  heat  f ux  is  proportional  to  (in)  . The  phase  is 
negative  for  all  frequencies  in  spite  of  the  low  frequency  amplitude  rise. 
The  ratio  of  the  perturbation-to-steady  local  heat  flux  has  been  plotted 
in  Figure  6(c)  against  a frequency  axis  scaled  by 

<1  .+■  (23) 

I + 2/Pr 


Except  for  some  discrepancies  associated  with  the  various  low  frequency 
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Figure  6(a).  The  local  perturbation  heat  flux  response  to  speed  modulations 
normalized  by  the  steady  heat  flux;  logarithmic  plot. 
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Figure  6(b).  The  local  perturbation  heat  flux  response  to  speed  modulations, 
normalized  by  the  steady  heat  flux;  linear  plot  detailing  the 
heat  f I ux  max i ma . 
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The  local  perturbation  heat  flux  response  to  speed  modulations 
normalized  by  the  steady  heat  flux;  a frequency  scaled  version 
of  (a ) . 
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peaks,  both  the  magnitude  and  phase  of  the  heat  flux  for  the  various 
Prandtl  numbers  satisfactorily  collapse  into  one  curve.  The  3 dB  band- 
width of  the  local  perturbation  heat  flux  is 


Q 


3dB 


1.45  { 


I + 2/Pr 
(I  + /Pr)2 


} 


heat  flux, 
(Schl ichti 


In  practice  one  senses  only  the  total  heat  flux.  The  total 
per  unit  width,  from  a plate  wetted  on  both  sides  is 


ng  1955,  pp.  268-269) 

= 4k  rJ  [f  ’ ’ ( 0)  ] 
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where  Rg  = U^L/v.  Similarly  the  total  perturbation  heat  flux  is 
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where  £ = wL/U^.  The  ratio  of  the  perturbat ion-to-steady  total  heat  flux 


2c[f"(0)]  V 


T.  ' (O.fi.PrJd^ 


is  plotted  in  Figure  7(a).  Asymptotically,  the  total  heat  flux  is  propor- 
tional to  (i£)  2,  in  contrast  to  the  local  heat  flux  ratio  which  is 
asymptotically  proportional  to  (ift)  ' . The  peaks  are  less  pronounced 
than  the  local  heat  flux  peaks.  The  Prandtl  number  scaling  of  equation 
(23)  satisfactorily  collapses  the  curves  for  various  Prandtl  numbers  into 


■f  1.0 


Figure  7(a).  The  perturbat ion-to-steady  total  heat  flux  ratio  for  speed 
modulations;  unsealed  frequency  axis. 
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one  curve  (Figure  7(b)).  The  3 dB  bandwidth  of  the  total  perturbation 
heat  flux,  the  anemometric  response  bandwidth,  is 
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3 dB 
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3b.  The  Temperature  Response 

Constant  +emperature  sensors  can  in  principle  be  used  to  measure 
temperature;  however,  in  practice  this  is  not  usually  the  case.  The 
effect  of  temperature  oscillations  of  the  ambient  fluid  is  in  practice  a 
contamination  of  the  anemometric  response.  The  transfer  function  of  the 
temperature  response  is  therefore  useful  for  estimating  the  contamination 
of  a measured  velocity  spectrum. 

The  simplified  version  of  the  temperature  response  equation  (21), 


name  I y 

e"  + ^ fPre'  - iqpre  = o 
0 = 0 § n ~ 0 


9=1  n ( 24 ) 

has  been  solved  by  the  same  routine  used  to  solve  the  momentum  equation  (12). 
The  numerical  grid  extended  to  n=IO  and  the  derived  temperature  field 
satisfies  equation  (24)  to  within  I part  in  I 06  or  better. 

The  local  heat  flux  response  to  variations  of  the  ambient  fluid 
temperature,  normalized  by  its  quasi-steady  value  is 

5T  e'(0,n,Pr) 

[f"(0)]Pr 

and  is  plotted  in  Figure  8(a),  for  the  Prandtl  number  range  0.7  to  15. 

In  contrast  to  the  anemometer  response,  the  local  heat  flux  response  to 
temperature  decreases  rapidly  with  increasing  frequency.  The  phase  shift 


1 
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and  the  rate  of  decrease  of  magnitude  (in  decibels  per  decade)  do  not 

tend  towards  an  asymptotic  limit.  The  curve  for  the  magnitude  of  the 

local  heat  flux  for  various  Prandtl  numbers  nearly  collapses  into  one 

curve  when  the  frequency  axis  is  scaled  by  Pr^'^^,  Figure  8(b).  The 

phase,  however,  does  not  have  a Prandtl  number  exponent  scaling.  The 

phase  shift  is  large,  even  at  very  low  frequencies.  Consequently  there 

are  regions  on  the  plate  where  the  heat  flux  direction  is  opposite  to 

that  at  the  leading  edge,  where  the  phase  shift  is  zero.. 

The  total  heat  flux  response  to  variations  of  the  ambient 

fluid's  temperature,  normalized  by  its  quasi-steady  value,  is 

x 

F J_ 

[f"(0)]Pr  5* 

0 

and  is  plotted  in  Figure  9(a).  The  character  of  the  curves  changes  around 
5=1.  In  this  frequency  region  the  local  heat  flux  phase  delay  near  the 
trailing  edge  of  the  plate  exceeds  270°.  This  local  "in  phase"  flux 
moderates  the  rate  of  flux  decrease  with  increasing  frequency.  For 
frequencies  above  C=l  the  relative  local  heat  flux  magnitude  is,  however, 
quite  small  and  consequently  the  total  heat  flux  is  pronortional  to 
( i£)-^  for  £>> I . 

A satisfactory  frequency  scaling,  for  both  the  magnitude  and 
0 38 

phase,  is  given  by  Pr  , Figure  9(b).  There  is,  however,  some  spreading 

0 38 

of  the  curves  near  5Pr  ’ =1  where  the  character  of  the  individual  curves 

depends  strongly  on  the  Prandtl  number. 

The  3 dB  bandwidths  for  the  response  to  temperature  oscillations 
in  the  main  stream  is 

— .zr~rr ....  _ ~ j 


e' (o,n,Pr)d^n 
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Figure  9(b).  The  total  heat  flux  response  to  a temperature  modulation  of 
the  ambient  fluid,  normalized  by  its  quasi~steady  value; 
scaled  frequency  axis. 


wLb  0.380 


This  makes  the  radian  frequency  temperature  band-width  0.56  U /L  at  Pr=l 

o 


and  0.25  IWL  at  Pr=8.  In  contrast,  the  corresponding  values  for  the 
anemometric  frequency  response  are  2.8  U^/L  and  1.7  Uq/L,  larger  by  a 
factor  of  5 and  6.8  respectively. 
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4.  Cone  I us  ions 

The  unsteady  viscous  stress  computed  numerically  from  a simplified 
unsteady  momentum  equation  that  is  not  restricted  by  frequency  compares 
favourably  with  the  frequency  restricted  results  of  Cheng  and  Elliott 
( 1957)  and  Lighthi I I (1954). 

When  the  boundary  layer  is  unsteady  (wx/Uq>> I ) the  equation 
governing  the  heat  flux  from  the  plate  in  response  to  speed  oscillations  in 
the  main  streams  is  significantly  different  from  the  equation  governing 
the  response  to  temperature  oscillations  in  the  main  stream.  Both  equations 
are  also  different  from  the  equation  governing  the  response  to  oscillations 
of  the  plate's  temperature.  The  un^feady  heating  source  that  generates 
the  anemometric  response  is  intrinsically  distributed  throughout  the 
boundary  layer  and  has  a profile  that  depends  strongly  on  the  Prandtl 
number  and  the  frequency  of  the  speed  oscillation.  The  heating  source 
that  generates  the  temperature  response,  on  the  other  hand,  is  constant 
for  all  frequencies  and  is  at  the  boundary  y=c°  for  external  temperature 
oscillations  and  at  y=0  for  internal  temperature  osc i I I at i ons . 

The  anemometric  heat  flux  response  increases  in  magnitude  with 
increasing  frequency  (from  its  quasi-steady  value  of  0.5)  to  a maximum 
value.  The  heat  flux  then  decreases  with  increasing  frequency  to  an 
asymptotic  value  proportional  to  (in)  ' for  the  local  heat  flux  and 

_ I 

(i£)~5  for  the  total  heat  flux  from  the  plate.  The  location  and  the 
value  of  the  heat  flux  maxima  depend  on  the  Prandtl  number.  The  largest 
maximum  occurs  at  a Prandtl  number  of  2.  The  phase  of  the  heat  flux  is 
always  negative. 
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The  local  heat  flux  in  response  to  external  temperature 

oscillations  decreases  with  increasing  frequency  and  has  neither  an 

asymptotic  logarithmic  rate  of  decrease  nor  an  asymptotic  phase  shift 

at  high  frequencies.  The  total  heat  flux  in  response  to  temperature 

oscillations  decreases  rapidly  from  its  quasi-steady  value  with  in- 

0 38 

creasing  frequency  up  to  the  frequency  £=1-4  Pr  and  then  becomes 

0 38  ~ - 

proportional  to  (£Pr  ' ) 2 for  larger  frequencies. 

The  -3  dB  bandwidth  is  3. 7(  I +2/Pr )/( I +/Pr  )■’  for  the  total 

_ Q 

heat  flux  in  response  to  speed  oscillations  and  is  0.56Pr  ' for  the 

total  heat  flux  in  response  to  temperature  oscillations  in  the  main 
stream. 


Frequency  response  calibrations  made  by  internal  heating  methods 
cannot  correctly  measure  the  response  to  speed  or  temDerature  oscillations 
in  the  main  stream,  for  constant  temperature  anemometers  or  thermometers, 
when  the  boundary  layer  is  not  auas i -steady . 
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Abstract 

The  steady  flow  calibration  data  of  two  platinum  thin  film 
turbulence  probes  show  that  their  thermal  conductance  Q/AT  is  not  pro- 
portional to  U 2 as  expected  for  a constant  temperature  conical  surface. 

A thermal  model  of  the  probes  that  agrees  favourably  with  the  calibration 
data  indicates  that  this  non-proport iona I i ty  is  due  to  the  thermal  effect 
of  the  glass  coating  over  the  film  and  glass  substrate  below  it.  The 
unsteady  heat  flux  from  the  film  in  response  to  variations  in  the 
ambient  fluid's  speed  increases  with  increasing  frequency  above  100 
Hertz.  It  is  shown  that  the  unsteady  heat  flux  is  proportional  to 
the  unsteady  viscous  surface  stress  over  the  film  which  increases  in 
magnitude  with  increasing  value  of  the  frequency  parameter  wL/311^. 

The  ratio  of  the  temperature-to-ve I oc i ty  sensitivity  is  highest  at 
zero  frequency. 


I ntroduct  i on 


The  steady  and  unsteady  calibration  of  two  (similar)  conical 
platinum  thin  film  turbulence  sensors  used  for  oceanic  measurements  are 
discussed  in  this  paper.  Comparable  probes  have  been  used  by  Grant, 

Stewart  and  ‘loilliet  (1962)  to  derive  the  Kolmogorov  spectrum  fron 
direct  turbulence  measurements  in  a tidal  channel,  and  by  Grant,  Moil  Met 
and  Vogel  (1968)  to  measure  turbulence  in  and  above  the  thermocMne. 

More  recently,  the  same  probes  discussed  here  have  been  mounted  on  a 
towed  bod/  for  open  ocean  measurements  (Nasmyth,  1970)  and  have  been 
mounted  on  a small  submersible  for  coastal  and  estuarine  studies 
(Gargett,  personal  communication). 

The  most  outstanding  feature  of  the  conical  hot  film  probe 
is  its  high  frequency  response.  The  frequency  response  is  calibrated 
from  10  to  1000  Hertz  by  vibrating  the  probe  in  the  direction  of  the 
steady  flow  in  a low  turbulence  water*  tunnel.  All  of  the  calibration  data 
presented  here  have  been  taken  by  the  Ocean  Mixing  Group  of  the  Institute 
of  Ocean  Science  at  Patricia  Bay,  British  Columbia,  Canada.  In  calibra- 
tion, the  modulation  in  velocity,  as  seen  by  the  probe,  arises  from  the 
oscillation  of  the  probe  parallel  to  a steady  stream.  However,  during 
field  measurements  the  modulation  in  velocity  stems  from  the  oscillation 
of  the  ambient  fluid  speed  parallel  to  the  steady  motion  of  the  probe. 

This  flow  difference  is  not  important;  the  behaviour  of  the  boundary 
layer  is  the  same  in  both  cases  when  the  fluid  is  incompressible 
(Lighthill,  1954).  Hence,  calibration  by  probe  vibration  correctly  simulate 
the  probe's  frequency  response  to  ambient  fluid  speed  oscillations. 
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It  is  more  common  to  measure  the  frequency  response  of  heated 
sensors  by  electronically  modulating  its  internal  temperature  in  the 
presence  of  steady  streaming  (Freymuth,  1977,  1978).  However,  this  method 
does  not  measure  a probe's  frequency  response  when  the  boundary  layer  is 
not  quasi-steady  (Lueck,  1979).  The  unsteady  calibration  data  presented 
indicates  that  the  boundary  layer  is  not  quasi-steady  at  frequencies 
above  approximately  |Q0  Hertz. 

A reproduction  of  the  probes  is  in  Figure  I.  The  basis  of  the 
probe  is  a blunt  glass  cone  having  a half-angle  of  15  degrees.  The 
sensing  element  is  a narrow  thin  platinum  film  which  is  deposited  around 
the  entire  circumference  at  a position  'downstream'  from  the  cone's  apex. 
Electrical  contact  with  the  film  is  made  on  opposite  sides  of  the  cone 
via  two  platinum  sponges  in  the  glass  substrate.  That  portion  of  the 
film  that  overlays  the  sponges  (hashed  area  in  Figure  I)  is  not  considered 
as  part  of  the  film's  sensing  surface.  A thin  quartz  layer  of  unknown 
thickness  covers  the  entire  nrobe  to  provide  electrical  insulation  from 
the  fluid  being  probed.  The  probes  were  made  under  special  contract 
at  the  University  of  British  Columbia.  The  probes  are  operated  in  a 
constant  temperature  mode  by  an  AC  servo-amp  I i f i er  centered  at  12.5  kHz. 

A detailed  account  of  the  probe's  electronics  and  its  operation  has  been 
given  by  Evans  (1962). 

The  steady  flow  calibration  curves  for  the  two  probes  ( V—  3 1 and 
V-38)  are  given  in  Figure  2.  Plotted  is  the  probe's  thermal  conductance 
Q/AT  against  the  square  root  of  the  Reynolds  number  where  0 is  the 
electric  power  dissipated  in  the  film  and  AT  is  the  temperature  difference 


Hot  Film  T-  AT 


Film  Surface  Area  s 0.92  (2x)Ld  sin <t> 

■ l.5*IO_3cm2 


V -38 


Film  Surface  Area  « 0.94(2T)Ldsin£ 

A sketch  of  probes  V-31  and  V-38,  their  relevant  dimensions  and  the 
coordinate  system  employed  in  the  analyses. 
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The  probes'  thermal  conductance  0/ AT  versur 
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between  the  film  and  the  ambient  water.  The  Reynolds  number  is  defined 
U L 

as  — — where  Uq  is  the  steady  ambient  fluid  speed  relative  to  the  probe, 

L the  distance  from  the  hypothetical  apex  to  the  start  of  the  film  (see 
Figure  I)  and  v is  the  kinematic  viscosity  of  the  ambient  fluid. 

The  heat  flux  from  an  ideal  constant  temperature  axi-symmetric 
surface  is  proportional  to  the  square  root  of  the  Reynolds  number  (White, 
1974,  p.  282)  if  the  usual  boundary  layer  approximations  are  valid. 

Bellhouse  and  Schultz  (1967)  reported  a linear  relationship  between  the 
heat  flux  and  the  square  root  of  the  ambient  fluid's  speed  for  a two- 
dimensional  sensor  in  air.  However,  the  data  in  Figure  2 show  curvature 
at  large  Reynolds  numbers  and  a systematic  deviation  away  from  a straight 
line  for  Reynolds  numbers  below  500.  This  deviation  at  lower  Reynolds 
numbers  cannot  be  attributed  to  a fundamental  change  in  the  Nusse't 
number  - Reynolds  number  relationship  of  the  flow.  Although  changes  in 
the  Nusselt  number  - Reynolds  number  relationship  have  been  reported 
(Hinze,  1975),  they  appear  to  be  restricted  to  Reynolds  numbers  below 
0.5.  In  addition,  Kramers  (1946)  reported  for  spheres  a linear  relation- 
ship between  the  Nusselt  number  and  the  square  root  of  the  Reynolds 
number  to  Reynolds  numbers  as  low  as  50.  Tan-atichat,  Nagib  and  Pluister 
(1973)  had  similar  difficulties  with  fitting  the  heat  flux  to  the  square 
root  of  the  speed,  for  a coated  cylindrical  sensor  in  water.  They  found 
that  plotting  the  heat  flux  against  U™  where  m = 0.225  gave  a moderately 
successful  fit  to  the  data;  however,  no  physical  explanation  was  given 
for  their  choice  of  m.  An  explanation  for  the  observed  steady  flow  response 


.k. 


is  presented  in  Section  2. 
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A plot  of  the  unsteady  response  of  probe  V- 38  is  shown  in  Figure 
3.  Plotted  is  the  probe's  unsteady  vol tage-to-speed  amplitude  ratio  in 
arbitrary  units  against  the  frequency,  in  Hertz,  of  the  speed  oscillation. 
The  lines  joining  the  data  points  are  a casual  fairing  of  the  data.  The 
unsteady  response,  at  a given  steady  speed,  decreases  slightly  with 
increasing  frequency  between  10  and  50  Hertz  and  then  increases  sharply 
with  increasing  frequency  up  to  1000  Hz.  The  character  of  the  frequency 
response  depends  strongly  on  the  ambient  fluid's  speed.  A more  significant 
response  decrease  at  low  frequencies  has  been  reported  by  Bellhouse  and 
Schultz  (1967)  for  a two-dimensional  probe  in  air.  These  authors 
attributed  the  response  decrease  to  an  attenuation  of  heat  waves  through 
the  probe's  substrate  and  predicted  the  effect  to  be  much  smaller  in 
water  than  in  air.  A small  difference  between  the  quasi-steady  response 
of  a hot  wire  anemometer  calculated  by  differentiating  the  steady  response 
and  the  low  frequency  response  observed  by  vibrating  the  probe  has  also 
been  reported  by  Morrison,  Perry  and  Samuel  (1972).  Their  calibrations, 
however,  were  restricted  to  one  frequency. 

An  analysis  of  the  bridge  circuitry  by  Evans  (1963)  indicates 
that  the  increase  in  response  with  increasing  frequency  above  50  Hz  cannot 
be  attributed  to  the  probe's  electronics  and  must  therefore  represent  a 
real  hydrodynamic  effect.  Evans  (1963)  showed  that  the  heat  flux  at  high 
frequencies  may  be  connected  to  the  unsteady  surface  shear.  The  nature 
of  the  probe's  unsteady  response  is  examined  in  Section  3. 


FREQUENCE  OF  SPEED  MODULATION  HERTZ 
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2.  Steady  Response 

The  steady  thermal  network  connecting  the  constant  temperature 
film  to  the  amb i ent  f I u i d w i I I be  mode  lied  by 

1)  a resistance  due  to  the  glass  coating  over  the  film  which  is  in  series 
w i th 

2)  the  thermal  resistance  of  the  fluid  boundary  layer  over  the  film  and  by 

3)  a shunt  resistance  through  the  glass  substrate  which  is  parallel  to  the 
coating  and  the  boundary  layer  resistance. 

The  electrical  analogue  of  the  proposed  thermal  network  is  shown  j 

in  Figure  4 where  the  subscripts  f,  c and  s refer  to  the  fluid,  the  coating  I 

and  the  substrate  respectively.  Because  the  heat  flux  through  the  substrate 
Qs  eventually  enters  the  ambient  fluid,  it  is  to  be  expected  that  this  heat 
flux  also  encounters  some  thermal  resistance  from  the  boundary  layer 
denoted  by  in  Figure  4.  The  slight  decrease  in  the  freauency  response 
between  10  and  50  Hertz  is  evidence  for  a small  boundary  layer  effect  in 
the  substrate  heat  path.  However,  the  bulk  of  the  thermal  resistance  of 
the  substrate  must  come  from  the  region  near  the  film  where  the  flux  lines 
are  relatively  concentrated.  The  thermal  resistance  R^'  is  determined 
by  the  surface  area  through  which  the  substrate's  heat  flux  Qs  passes 

• J 

upon  entering  the  fluid,  by  the  thermal  conductivity  of  the  fluid  and  by  the 
efficiency  of  forced  convective  heat  transfer  in  the  boundary  layer.  Because 
this  surface  area  is  comparable  to  or  greater  than  the  surface  area  of  the 
film  R^'sRf.  Moreover,  the  conductivity  of  water  at  20°C  (,v5.9xl0"3  watts/ 

°C-cm)  is  comparable  to  the  conductivity  of  glass  ('vl.QxlO-2  watts/°C-cm) 
whereas  the  conductivity  of  air  is  much  smaller  (only  ,''2.5xl0-4  watts/  C-cm). 


Subsequent  analyses  show  that  for  the  probes  in  question  R^<<Rs  anc*  hence 
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R,’«R  in  water.  In  air,  however,  it  is  possible  that  R/'vR  . In  the 
f S 

remaining  work  Rf'  is  ignored  relative  to  Rs>  the  substrate's  resistance. 
The  probe's  coating  resistance  is  approximately 
D 


R = , c 
c X S 
c 


(I  ) 


where  D is  the  coating  thickness,  Xc  is  its  thermal  conductivity  and 
S is  the  surface  area  of  the  film.  The  resistance  of  the  boundary  layer 
can  be  written  as 


R,  = — 

f /IT 


(2) 


if  it  is  dominated  by  forced  convection,  where  Rg  is  the  Reynolds  number 
and  a,  which  will  be  derived  later,  is  a flow  independent  constant  that 
depends  on  the  properties  of  the  fluid  and  the  probe's  geometry.  The 
appropriateness  of  equation  2 has  been  demonstrated  by  White  (1974,  p.  282). 

The  equivalent  impedance  of  the  circuit  in  Figure  4 with  R^'=0  is 

B 


Z = ^ = C + 


A + /R~ 
e 


where 


C = 


R R 
c s 

R + R ’ 
c s 


(3) 


(4a) 
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number  equation  3 gives  a thermal  resistance  equal  to  the  parallel 
equivalent  of  the  coating  and  substrate  resistance  and  at  very  small 
Reynolds  numbers  a resistance  equal  to  the  substrate  resistance.  Thermal 
convection,  which  must  be  present  when  Rg< I , is  not  considered. 

The  available  steady  flow  calibration  data  for  probes  V— 3 1 and 
V-38  has  been  fitted  to  equation  3 by  a least  square  error  routine 
to  determine  the  constants  A,  B and  C.  All  measurements  were  made 
at  an  ambient  water  temperature  of  20°C.  The  data  are  plotted  in 
Figure  5.  The  solid  line  is  the  thermal  resistance  Z=AT/Q  determined 
by  the  least  square  routine  and  the  symbols  represent  the  observed 
resistances.  The  quality  of  the  fit  in  terms  of  the  percent  difference 
between  the  fit  and  the  observed  values  is  plotted  near  the  top  of  the 
graphs.  The  fitted  impedance  of  the  probes  are 


f ■ ^.4  * 


3220 


7.77  + /if 
e 


C/watt 


(5a) 


for  V— 3 1 and 

41  = 12.5  + °C/watt  (5b) 

V 12.3  + v^T 

e 

for  V-38.  Table  I summarizes  the  experimentally  determined  values  of 

R , R and  a.  The  data  of  zero  flow  were  not  used  to  calculate  the  least 
s c 

square  fit.  Some  care  was  taken  to  attempt  to  minimize  thermal  convection 
at  zero  flow  speeds.  Al I of  the  data  for  V— 3 1 agree  to  better  than  4% 
and  all  except  three  of  the  lower  speed  points  agree  to  better  than  2% 
with  the  fitted  curve.  All  of  the  data  for  V-38  agrees  to  better  than 
0.8$  with  the  fitted  curve.  The  two  data  sets  for  V-38  were  taken  one 


4 
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5(a)  The  probes’  steady  thermal  resistance  AT/Q  versus  Rg2.  The  solid 
lines  represent  the  least  square  fit  of  the  thermal  model  to  the 
data.  The  quality  of  the  fit  is  indicated  near  the  top  of  the 
figure;  for  V-31 . 
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7.  ERROR 
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V-3 1 


A = 7.77 
B = 3220  °C/watt 
C = 25.4  °C/wa+t 
S = (0. 92)  2irLds  i n4> 


R = 440  °C/watt 
s 

R = 27  °C/watt 

c 

a = 3630  °C/watt 
1 . 50  x 1 0-3  cm2 


V-38 


A = 12.3 

R 

s 

B = 4410  °C/watt 

Rc 

C = 12.5  °C/watt 

a 

S = (0. 94)2irLds  i nt  = 

1.20  x 1 0-3  cm2 

370  °C/watt 
13  °C/wa+t 
4720  °C/watt 


Tab  I e B I 

Thermal  properties  of  the  hot  film  probes  derived  from  fitting  the 
thermal  model  to  the  steady  calibration  data. 
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month  apart  and  show  a systematic  difference  of  approximately  0.8/5.  It 
appears  that  the  difficulty  with  fitting  the  probe's  heat  flux  to  a 
proportionate  relationship  with  is  due  to  the  thermal  effect  of  the 
coating  and  the  substrate. 

The  coating  thickness  over  the  probes,  using  equation  (I)  and  the 
surface  areas  listed  in  Table  I,  is  4.1  x 1 0-4  cm  for  V— 3 1 and 
1.6  x 1 0 4 cm  for  V-38.  The  coating  thickness  on  these  probes  is  other- 
wise unknown  but  two  commercial  manufacturers  of  similar  probes  coat 
theirs  to  a nominal  thickness  of  2 x I0-4  cm.  The  substrate  resistance 
should  be  given  by 
D' 


R = . c 
s X S 
s 


(6) 


where  D'  is  a character i st  ic  depth  of  thermal  penetration  and  should  be 
comparable  to  the  width  of  the  film.  Using  equation  (6)  gives 
D'  = 67  x I CT4  cm  for  V— 3 1 which  is  almost  exactly  equal  to  the  width 
of  the  film  and  D’  = 46  x 10  4 cm  for  V-38  which  is  20$  less  than  the 
film's  w i dth. 

2(a).  Boundary  Layer  Resistance 

The  hydrodynamic  equations  for  the  probe's  boundary  layer,  using 
the  coordinates  shown  in  Figure  I,  are  (White,  1974,  p.  340) 


f-  (Ru  ) + R v =0 
3x  o 3y  o 


3u  3u 

o , o 

U + v T = V 

o 3x  o 3y 


32u  dV 

— ; r + v 

3y2  ° dX 


(7) 


(8) 


u(x,y)  = v (x ,y ) = 0,  u(x,®) 

0 0 0 


= V (x) 
0 


1 
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where  VQ  is  the  steady  potential  flow  at  the  probe's  surface,  uq  and 

v are  the  steady  boundary  layer  velocities  and  R = R(x)  is  the  radius 
0 

of  the  probe.  The  potential  flow  over  an  ideal  axi symmetric  cone  is 
(White,  1974,  p.  341 ) 


V (x)  = Ex 

o 


(9) 


where  E is  a constant  and  n=0.035  for  a cone  of  half-angle  15  . Equation 
8 reduces  to  the  ordinary  differential  equation  (White,  1974,  p.  242) 

2n  (l-f'2)  = 0 


, » 1 1 ii 

f + ff  + 


3+n 


f ' (0)  = f (0)  = 0 


f (») 


(10) 


where 


u (x,y)  = V (x)  f (q)  = ExP  f'(q) 
o ' o 


(II) 


, . /2 c n- 1 r n+3  , . n- 1 , * , 

/o(x’y)  = _ /3+n  Ex  {~  f + ~T  nf  } 


(12) 


A 


/3+n  c n- I 

n = y /^r—  Ex 


2v 


and  prime  denotes  differentiation  with  respect  to  q.  Equation  10  is 
identical  to  the  equation  for  two-dimensional  flow  over  a wedge  of  half- 
angle 2.1  degrees.  Solutions  to  equation  (10)  for  various  values  of  q 
can  be  found  in  White  (1974,  p.  276)  and  Schlichting  (1955,  p.  129). 

At  the  blunt  tip  of  the  cone  and  at  the  rear  of  the  conical  section 
the  potential  flow  at  the  surface  cannot  be  represented  by  equation  9. 

Over  most  of  the  conical  section  of  the  probe  equation  9 should  be  accurate, 


and  because  n=0.035<<l  the  potential  flow  Vq  should  not  differ  significant- 


ly from  the  free-stream  speed  U . However,  a detailed  analysis  of  the 


potential  flow  over  probes  V— 3 1 and  V-38  has  not  been  made. 
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The  equation  governing  the  boundary  layer  temperature  TQ(x,y),  measured 

with  respect  to  the  ambient  fluid's  steady  temperature,  is  (Lighthill,  1954) 

30  A 30  320  . - ... 

u -r-  + v -r K — - =0  (14) 

o 3x  o 3y  „ 2 

3y 

where  0 = ATq  - TQ(x,y),  ATq  is  the  temperature  difference  between  the 
probe's  wetted  surface  and  the  ambient  fluid,  and  k is  the  thermal 
diffusivity  of  the  fluid.  Assuming  negligible  heating  upstream  of  the 
film  and  a uniform  temperature  ATq  at  the  wetted  surface  over  the  film, 
the  boundary  conditions  on  equation  14  are 

0 = 0 at  y = 0 L$x$L+d 

0 = AT  for  a I I x<L  ■ (15) 

o 

0 = AT  as  y-*»  for  all  x. 
o 7 

Because  the  momentum  boundary  layer  starts  at  the  tip  of  the  cone  whereas 
the  thermal  boundary  layer  starts  at  the  film  (x  = L)  the  similarity 
variable  n is  not  appropriate  for  the  temperature  equation  (14). 

Following  Kays  (1966)  let  the  velocity  profile  be  approximated  by 


the  cubic  parabola 


u (x,y) 


i “0  ’r  3 v I v 

f (n)  = - = 4 < 4 ) - ~ ( 4 


V (x) 
o 


^ ( JL  ) _ ( JL  ) 

2 6 2 6 


where  6 is  an  artificial  viscous  boundary  layer  thickness.  The  choice 


6 = 3.04 


/ W- 

(3+n)  Ex 


makes  the  surface  velocity  gradient  calculated  from  equation  (16)  agree 
with  the  value  tabulated  in  White  (1974,  p.  270).  Assuming  a similar 


profile  for  temperature  we  can  write 
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6 _ 3 , y_  x I , y 

at-  = j ( h } ■ 2 ( h ) (l8) 

0 I T 

where  6^  is  the  artificial  thermal  boundary  layer  thickness.  The  thermal- 
to-viscous  boundary  layer  thickness  ratio  r = 6T/6  is  (Kays,  1966,  p.  218) 


( - )C  } 
X 


where 


f ( ^ ) “ 3/4 


and  Pr  = v/k  is  the  Prandtl  number. 


The  steady  heat  flux  from  the  probe's  surface  can  be  calculated 
using  equations  17,  18  and  19,  and  is 


/ \ •>  59 

V*’  ’ xf  57 


, 4T  fr—  . A.Pr  /U  > 

2 o__  /3+n  J_  f / _o_ 

2 3.04/  2 x r|_(i)Ci|/3  v 

y=0  ‘ x ^ 


where  is  the  thermal  conductivity  of  the  fluid.  The  variation  of 
the  heat  flux  over  the  wetted  surface  is  dominated  by  the  factor 

L c 1/3 
( I - (-)c  ) 
x 

because  the  film  is  relatively  narrow  (d/L  <<  I).  Accordingly,  x can 
be  set  equal  to  L where  it  occurs  outside  of  the  above  factor  and  this 
factor  can  be  set  equal  to  (c8  / L)1^  where  6 = x - L <<  L.  The  total 
heat  flux  into  the  fluid  above  the  film  is  then  obtained  by  integrating 
equation  21  over  the  surface  area  of  the  film; 


TU 

Q0  = S j QQ(x)dx  = ( ) 2itLs i n<J>  QQ(x)dx  (22) 

L L 

_ 9 , 0.92  . 2TrLsin<|)  ,3+n. 1/2  _4, 1/3 fd  2/3  { 1/3 

" 4 ( 0.94  > 3T04  (~5  3-Hn  ^ We  Pr  * 


(x)dx  (22) 


r 


The  thermal  resistance  coefficient,  equation  2,  is 


9o  9 .0.92  2nLs  i n<f>  ,3+n  l/2,4_)l/3(d,2/3  p 1/3 
0 1 = ~ TT  = 4 (0.94)  ^04  ( 2 5 W Vr  ' 


at  r 

o e 


(23a) 


which  for  both  probes  equals  (to  within  \%) 


d 2/3  1/3 

.53  Xf  L (^)  Pr 


(23b) 


Using  the  probes'  dimensions  (figure  I),  = 6.3  * I0-3  watts/ 

°C-cm  (appropriate  for  water  at  40°C),  and  a Prandtl  number  of  7.1 
(water  at  20°C)  in  equation  23  predicts  values  for  a of  2880°C/watt 
for  V— 3 1 and  3320°C/watt  for  V-38.  The  experimentally  derived  values 
for  a are  3630  and  4720  °C/watt.  The  boundary  layer  heat  flux  infered 
from  fitting  the  thermal  model  to  the  calibration  data  is  smaller  than 
the  calculated  heat  flux  by  2 1 % tor  V—  3 1 and  30?  for  V-38.  Some  of 
this  discrepancy  is  due  to  a variable  Prandtl  number  in  the  thermal 
bounday  layer.  At  40°C  the  Prandtl  number  is  4.3,  for  instance.  In 
addition  the  origin  of  the  momentum  boundary  layer  has  been  assigned 
to  the  hypothetical  apex  (x=0)  rather  than  the  real  (blunted)  apex. 

Placing  the  origin  of  the  momentum  boundary  layer  at  the  real  apex 
would  accentuate  the  discrepancy  whereas  decreasing  the  Prandtl  number 
to  account  for  the  warming  of  the  thermal  boundary  layer  would  diminish 
the  discrepancy. 

Some  implications  of  the  proposed  thermal  model  can  now  be  examined. 
Oceanic  measurements  with  V — 3 1 and  V-38  are  typically  taken  at  130  cm/sec. 
This  corresponds  to  a Reynolds  number  of  approximately  1800.  A typical 
value  for  the  boundary  layer  resistance  at  that  speed  is 
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'I 

Rt  = — * 85°C/watt. 

f 

e 

Takinq  R = 27  and  R = 440°C/watt  (Table  I)  one  calculates  the  temperature 
3 c s 

drop  across  the  coating,  AT^,  from 

1TC  ■ 0.24AT. 

f c 

Twenty-four  percent  of  the  f i Im-to-ambi ent  fluid  temperature  drop  occurs 
across  the  coating.  The  boundary  layer  resistance  coefficient,  a,  is 
considerably  larger  in  air  than  in  water.  Placing  the  appropriate  values 
of  fluid  conductivity,  A^,  and  Prandtl  number  into  equation  23  gives 

a . = I . 7 x I 05  °C/watt 

a i r 

A Reynolds  number  based  on  a more  representative  air  flow  speed  of  10 
meters  per  second  gives 

Rf  = 4500  °C/watt. 
a i r 

According  to  the  model  the  temperature  drop  across  the  coating  in  air 
then  i s 


AT 

c 


AT 


27 

27+4500 


6 x |0-3  AT. 


Only  0.6  percent  of  the  f i Im-to-ambi ent  fluid  temperature  drop  occurs 
across  the  coating  when  the  same  probes  are  operated  in  air.  The 
coating's  thickness  and  the  coating's  thermal  conductivity  are  obviously 
much  more  important  to  a probe's  heat  flux  response  in  water  than  in 
a i 


A fraction  of  the  heat  flux  from  the  film  goes  through  the  substrate 
and  consequently  contributes  nothing,  if  R^<<Rg,  to  the  probe’s  speed 
sensitivity.  Let  denote  the  substrate  flux  and  Qq  the  boundary  layer 
flux  (Figure  4).  The  flux  ratio  is 


R,  + R 
f c 


For  water  at  Uq  = 130  cm/sec,  Qs/Q0  is  approximately  0.25;  i.e.  20%  of 

the  total  flux  from  the  film  goes  into  the  substrate.  In  contrast, 

Qs/Qq  should  be  approximately  14  in  air,  indicating  that  more  than 

00%  of  the  total  flux  is  into  the  substrate.  However,  this  part  of 

the  model  is  not  completely  applicable  in  air  because  r|  may  be 

comparable  to  R . Nonetheless,  the  same  probes  operated  in  air 
s 

would  lose  a significant  amount  of  speed  sensitivity  because  of 
conduction  into  the  substrate.  A thin  film  turbulence  probe  suitable 
for  use  in  water  may  be  unsuitable  for  use  in  air. 

The  quasi-steady  sensitivity  to  both  speed  and  temperature 
fluctuations  in  the  ambient  fluid  are  now  readily  attained.  The 
quasi-steady  sensitivity  to  speed  is 

9Q  _ 92_  _^e  = ^e 

3U  3R  3U  U 3R 

o e o o e 

and  to  temperature  is 
9Q  _ 

a (AT)  AT  • 


The  ratio  of  the  quasi -steady  speed-to-temperature  sensitivity  is 
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9Q/9U 

0 

Re  90  , AT  . 

= ^ 9R  ( U~  ) 
e o 

9Q/9(AT) 

\ 


B R 


(A+R  * ) (AC+B+CR  * ) 
e e 


AT 

2U 


The  ideal  conical  probe  has  no  coating  resistance  and  an  infinite 

substrate  resistance,  i.e.  A = C = 0.  Thus  the  optimum  relative 

sensitivity  is  AT/2U  . The  real  probe  (V-31)  has  at  U =130  cm/sec 

o o 

a relative  sensitivity  of 

0.61  (§J_>  , 

0 

which  is  61$  of  the  optimum  sensitivity  at  130  cm/sec.  Using 
AT  = 20°C  gives 

9Q/9U 

9q7V(at)  = 0*047°c  per  cm/sec 

which  means  that  the  probe  cannot  distinguish  a quasi-steady  fluctuation 
of  0.047°C  from  a fluctuation  of  I cm/sec  at  130  cm/sec. 


3(a).  Viscous  Boundary  Layer 

If  the  ambient  fluid's  speed,  relative  to  the  probe,  has  an 
infinitesimal  sinusoidal  oscillation  about  its  steady  mean,  then  one 
can  write  the  potential  flow  at  the  probe's  surface  as 

V(x,t)  = V (x) ( l+ce'w^)  = Exn(l+ee'w^)  = U (l+ee'^)  (24) 

o o 

where  u>  is  the  angular  frequency  in  radians/sec  and  where  e,  small  compared 
to  unity,  is  the  normalized  speed  amplitude  of  the  probe  shaker  used  for 
response  calibrations.  The  boundary  layer  velocity  will  consist  of  a 
steady  mean  velocity  plus  a sinusoidal  perturbation,  i.e. 


(27) 


■ 3k  {2n  * 3,21  - 2(5^>  2 - f'9''  (29> 

with  boundary  conditions 

g(o,0)  = g'(o,C2)  = 0 

g ' (°°,0)  = I 

g(n,0)  = i (f  + nf ' ) . 

The  presence  of  n in  the  momentum  equation  (29)  results  from  the  steady 
pressure  gradient  along  the  cone's  surface.  Its  effect  relative  to  a 
zero  pressure  gradient  flow,  is  an  increase  in  the  fluid's  inertia  (the 
third  term  on  the  left-hand  side  in  (29))  and  an  increase  in  unsteady 
pressure  gradient  (the  first  term  on  the  right-hand  side).  The  steady 
pressure  gradient  along  a cone  of  half-angle  15°  is  evidently  quite 
small  because  n = 0.035.  Thus  n can  be  ignored  relative  to  3 and, 
because  0 s f ' s I , 2n  can  also  be  ignored  relative  to  3ift  when 


mmmmm 
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Q >>  2n/3  = 0.023.  Within  the  above  constraints,  the  perturbation 
momentum  equation  (29)  simplifies  to 

g'"  + fg"  - 2ifig'  + f"g  = -2  in  + § <1  In  (f ’’g-f’g* ) (30) 

g(0,0)  = g' (0,n)  = 0 

g'(o,,fi)  = | 

g(n,0)  = i(f+nf'). 

Except  for  the  definition  of  the  variables  n and  0 equation  (30)  is 
identical  to  the  corresponding  perturbation  momentum  equation  for 
two-dimensional  flow  over  a flat  plate  (Lueck,  1979).  The  solution  for 
the  unsteady  flow  over  a plate  also  applies  to  the  unsteady  flow  over 
a cone  of  half-angle  15°,  within  the  limits  cited. 
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The  artificial  concept  of  a finite  boundary  layer  thickness  will  be 

used  to  solve  T because  n cannot  be  a similarity  variable  for  T ^ , as 

was  the  case  for  T . 

o 

Let 

T j ( x , y , u ) = T j ( £ , ft ) 

where,  using  equations  (13),  (17)  and  (19) 

r - = = n.  . (33) 

^ 6 r6  3.04r 

T 

and 

0 < £ < I . 

The  steady  viscous  boundary  layer  thickness,  6,  is  essentially  constant 
over  the  film  because  the  aspect  ratio  d/L<< I . Hence,  x can  be  set 
equal  to  L except  where  it  appears  in  r.  The  differential  operators 
in  the  temperature  equation  (32)  then  become 

iL_  = — L iL_  = — ! /_°  <L_  (34a) 

3y  r6  3£  3.04r/  2vL  3£ 

2 2 3U  2 3 

3 ..  _ 1 9 = ! o i — ( 34b ) 

2 2^  v2  2\)L  2 

3y  (r6T  3 £ (3.04r)z  3? 
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Using  the  above  operators  along  with  equations  (II),  (12),  (18),  (19), 
(27)  and  (28)  transforms  the  terms  in  the  perturbation  temperature 


equation  (32)  to 
U 


(a)  7-—  x 3 i OT 

L 1 


(b)  - p-  x {( l+2r3Pr)£f ' - 


9T . 


4r 3Pr 


3?T, 


(c)  - ~ x 


2(3.04 ) 2r 2Pr  3£2 


(35a) 


(35b) 


(35c) 


U 3 ( AT  ) ( l-C  ) , , 

(d)  - x {(l+2r  Pr)£g  - 


3.04  <69+4!h»>-  <35d> 


For  the  two  probes  considered,  r is  always  less  than  0.17.  Thus  2r  Pr 
is  small  compared  to  unity.  Throughout  the  domain  0 < C < I the 
variable  n never  exceeds  0.51;  consequently,  the  steady  velocity 
profile  f'  is  essentially  proportional  to  n (Sch I icht inq , 1955, 
p.  107).  A linear  approximation  for  the  unsteady  velocity  profile  g 
is  reasonable  up  to  frequencies  0 not  much  greater  than  unity  (Lueck, 
1979).  Thus  for  n < i and  0 < I the  boundary  layer  velocity  profiles 
over  the  film  can  be  approximated  by 

f = i nf"(0)  = ^(3.04)2rVf " (0) 
f ’ = nf ' ’ (0)  = 3.04  r£f ’ ' (0) 

g = i ng"(0,0)  = i (3.04)2r2eV  (0,0) 
g'  = ng'  (0,0)  = 3.04  r£  g"  (0,0) 


Substituting  the  above  into  the  perturbation  temperature  equation  (35 
and  32 ) yields  finally 
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3 T1  (3  04) V * * 33 ] 2 2 

— y + ^ J-  * f (0)  — 1 - 2(3.04)  r Pr  i f>T 

z O oc,  1 


— C2(|-e)(&To)9"(0,Q)  +-0(r3Pr) 


(36) 


Tx  = 0 at  C = 0,  I 
2 2 

The  term  2(3.04)  r PriOT^  which  represents  the  thermal  inertia  of  the 
boundary  layer,  is  evidently  small  compared  to  the  other  terms  in 
equation  36  for  the  frequency  ranpe  considered  here  Os  I • Oroppinq 
the  third  term  in  equation  36  reduces  it  to  the  ordinary  differential 
equat i on 


dTl  . (3.04) 3 2 »'  n,  y _ ( 3. 04 ) 3 .2 

d^2  6 5 dC  “ 


5 (l-£  HAT  )g  (0,P) 

(37) 


T = 0 at  £ = 0, 


Equation  37  does  not  have  to  be  solved  explicitly.  The  perturbation 
temperature,  Tj,  is  everywhere  proportional  to  (AT  )g' ' (0,0) . Hence 
the  total  heat  flux  from  the  film  is 


0 = Qq  + q (0  )e  ' W+  = 0Q  + ee ' W+F(ATo )g ' ' (0,0) 


(38) 


where  F is  a frequency  independent  constant.  In  the  limit  0-*-0,  the 
perturbation  heat  flux  must  equal  the  quasi-steady  heat  flux,  i.e. 

9Q- 


4(0).  = M = 

e o 9U  ‘ o 


o 9u  i Q~ = F(y)g  (o'o) 


The  constant  F is  readily  determined  by  using 


» ! 3 i t 

g (0,0)  = | f (0). 
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The  total  heat  flux  into  the  fluid  above  the  film  then  is 


Q = Q [|  + | eiu+  3— 

0 J f (0) 


(39) 


or  using  (23b) 


h 2/3  1/3  1/2  c iut  n"  (0  il! 

Q = I .53X.L(f>  Pr  / R AT  [I  +|e  “ 9 — ] 
f L e 0 3 fM(0) 


(40) 


The  unsteady  response  of  the  two  thin  film  probes  is  plotted  in 

figure  6 for  the  steady  speeds  102,  129,  154  and  182  cm/sec.  Plotted 

1 

is  the  unsteady  thermal  conductance  q/ AT q divided  by  and  the 

fractional  amplitude  e against  the  frequency  parameter  wL/3Uq.  The 

unsteady  conductance  has  been  scaled  by  an  arbitrary  factor,  y,  to 

bring  the  low  frequency  values  close  to  0 dB.  (The  calibration  for 

the  velocity  monitor  on  the  shaker  is  unavailable.)  The  solid  line  in 

figure  6 is  2/3  g''(0,Q)  / f'(0)  which  accord i ng  to  equation  40  is 

1 

proportional  to  q/(e&T  U 2 ) . The  dashed  line  has  been  drawn  to  show  the 

0 0 

quality  of  the  fit  if  the  frequency  parameter  it>L'/3U^  had  been  used 

instead  of  ioL/3IJ  (figure  I).  The  data  for  both  V— 3 1 and  V-38  exhibit 
0 

a rise  with  increasing  frequency  characteristic  of  g''(0,R)  but  appear 
to  be  shifted  by  approximately  40$  to  the  right  along  the  frequency  axis 
(less  if  scaled  by  <dL?/3U  ).  This  shift  along  the  frequency  axis  may 
reflect  some  inadequacies  with  modelling  the  flow  over  a blunted  cone 
by  the  flow  over  an  ideal  "needle  sharp"  cone.  The  scatter  among  the 
data  points  for  R> I . 5 is  noticeable.  This  scatter  may  be  due  to  some 
spurious  motions  of  the  probe  support  at  high  frequencies  (Evans, 
personal  communication);  I . 5 corresponds  to  real  frequencies  of  500 
to  1000  Hertz.  Thermal  inertia  should  become  important  for  Q> I and 
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The  probes'  unsteady  thermal  conductance  scaled  by  U 3 and  the 
fractional  speed  amplitude  e versus  the  frequency  parameter 
0 = u>L/3U0.  The  solid  line  is  two-thirds  of  the  unsteady-to- 
steady  surface  stress  ratio  after  Lueck  (1979).  The  dashed  line 
indicates  the  quality  of  the  fit  if  the  frequency  had  been  scaled 
by  u)L'/3U0.  L'  is  shown  in  fiqure  I. 


cause  a "rol I -off"  in  the  frequency  response,  but  only  the  data  for 
Uq  = 1 02  cm/sec  shows  a definite  trend  to  roll-off.  Although  the  match 
between  the  data  and  the  theoretically  predicted  response  is  not 
complete  it  does  appear  that  the  unsteady  heat  flux  from  the  film  is 
strongly  influenced  by  the  unsteady  viscous  boundary  layer  and  that 
for  probes  V— 31  and  V-38  it  is  proportional  to  the  unsteady  surface 
shear  g'  ' ( 0,f) ) . 

The  peculiar  rise  of  the  response  with  increasing  frequency  is 
due  to  the  probes'  aspect  ratio  d/L«l  . Most  commercial  probes  have 
aspect  ratios  closer  to  unity  and  therefore  should  show  much  less  if 
any  rise,  at  high  frequencies  due  to  the  thermal  inertia  of  the  boundary 
layer.  Observations  by  Nowell  (1974)  of  commercial  constant  temperature 
hot  film  probes  show  a decrease  in  response  with  increasing  frequency 
above  apDrox i mate  I y 100  Hertz. 

The  sensitivity  of  these  probes  to  temperature  oscillations  in 
the  main  stream  causes  a contamination  of  the  measured  velocity  signal 
which  can  be  very  severe  (Gargett,  1978).  The  response  of  these 
probes  to  temperature  oscillations  is  difficult  to  calculate.  However, 
at  least  three  aspects  of  the  response  to  infinitesimal  external 
temperature  oscillations  can  be  deduced.  Firstly,  the  heat  flux  into 
the  boundary  layer  above  the  film  cannot  increase  with  increasing 
frequency  because  there  is  no  mechanism  to  accentuate  the  response  to 
temperature  oscillations  in  the  main  stream.  As  well,  calculations  for 
an  ideal  flat  plate  sensor  by  Lueck  (1979)  show  that  the  bandwidth  of 
the  response  to  temperature  oscillations  is  smaller  than  the  bandwidth 
for  the  response  to  speed  oscillations.  Secondly,  the  heat  flux  into 


n 
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the  substrate  induced  bv  external  temperature  oscillations  will  decrease 
with  increasing  frequency  because  of  the  substrate's  finite  diffusivity. 

The  effective  depth  of  penetration  into  the  substrate  is  approximately 
the  width  of  the  film  d = 6 x I0"3  cm,  hence  the  heat  flux  decreases 
when  ud2/Ks  * I,  i .e.  above  approximately  30  Hertz.  Thirdly,  the  alass 
ups i ream  from  the  film  acts  as  a passive  heat  sink  and  reduces  somewhat 
the  temperature  oscillations  of  the  fluid  close  to  the  surface  before 
this  fluid  reaches  the  area  above  the  film.  Thus  the  ratio  of  temoerature- 
to-velocitv  sensitivity  is  hiqhest  at  zero  freauency  and  decreases  with 
increasing  frequency.  An  upper  bound  to  the  contami nation  of  a measured 
velocity  spectrum  by  temperature  fluctuations  can  be  calculated  from 
the  ratio  of  ouasi-steady  temperature-to-ve I oc i tv  sensit;vity  and  an 
independent  estimate  of  the  temperature  spectrum. 


I 


4.  Cone  I us  ions 
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The  steady  calibration  data  presented  show  that  the  probes'  thermal 

I 

conductance  Q/AT  is  not  proportional  to  U*  as  expected  for  a constant 
temperature  conical  surface.  A thermal  model  of  the  probes  indicates  that 

I 

0/AT  is  neither  proportional  nor  linearly  related  to  U*  because  of  the 

thermal  effects  of  the  coating  over  the  film  and  the  substrate  below  it. 

The  model  predicts  that  a glass  coating  of  only  n.  | x 1 0-4  cm  can 

significantly  reduce  the  temperature  difference  between  the  probe's  wetted 

surface  and  the  a/nbient  fluid  and  that  the  heat  flux  into  the  substrate 

is  moderately  important  for  probes  oDerated  in  water.  For  probes 

operated  in  air,  on  the  other  hand,  the  model  predicts  the  opposite 

characteristics:  a relatively  small  effect  due  to  the  coating  and  a 

very  large  heat  flux  into  the  substrate.  The  coating  and  the  substrate 

reduce  the  quasi -steady  speed-to-temoerature  sensitivity  of  the  probes 

by  about  40$  from  its  optimum  value  of  AT/2U  . 

o 

The  unsteady  heat  flux  from  the  hot  film  probes  in  response  to  varia- 
tions of  the  ambient  fluid's  speed  is  stronily  influenced  by  the  unsteady 
viscous  boundary  layer.  For  the  probes  considered,  the  unsteady  heat 
flux  is  proportional  to  the  unsteady  viscous  surface  stress  over  the 
film.  The  unsteady  heat  flux  from  the  probes  examined  increases  in 
magnitude,  over  its  quasi-steady  value,  with  increasing  frequency.  The 
ratio  of  temperature- to-ve I oc i ty  sensitivity  is  highest  at  zero  frequency 
and  decreases  with  increasing  frequency.  An  upper  bound  to  the  contamina- 
tion of  a measured  velocity  spectrum  by  temperature  fluctuations  in  the  main 
stream  can  be  calculated  from  the  ratio  of  the  quasi-steady  temperature- to- 


velocity  sensitivity  and  an  independent  estimate  of  the  tempeature  spectrum. 


Bi b I ioqraphy 


Bellhouse,  B. J . and  D.L.  Schultz.  1967.  The  determination  of 

fluctuating  velocity  in  air  with  heated  thin  film  gauges. 
Journal  of  Fluid  Mechanics,  29g  289-295. 

Evans,  D.J.  1963.  An  instrument  for  the  measurement  of  ocean 

turbulence.  Technical  Memorandum  63-8,  Pacific  Naval  Lab- 
oratory, Defense  Research  Board,  Canada. 

Freymuth , P.  1977.  Freauency  response  and  electronic  testing  for 
constant-temperature  hot  wire  anemometers.  Journal  of 
Physics  E;  Scientific  Instruments,  10:  705-710. 

Freymuth,  P.  1978.  Theory  of  frequency  optimization  for  hot  film 

anemometers.  Journal  of  Physics  E;  Scientific  Instruments, 

I I : 177-179. 

Gargett,  A.E.  1978.  Microstructure  and  f i nestructure  in  an  upper 
ocean  frontal  regime.  Journal  of  Geophysical  Research, 

8 3 ( C I 0) : 5123-5134. 

Grant,  H.L.,  R.W.  Stewart  and  A.  Moil  Met.  1962.  Turbulence  spectra 

from  a tidal  channel.  Journal  of  Fluid  Mechanics,  J_2:  241-288. 

Grant,  H.L. , A.  Moil  Met  and  V.M.  Vogel.  1968.  Some  observations  of 
the  occurrence  of  turbulence  in  and  above  the  thermocline. 
Journal  of  Fluid  Mechanics,  34:  443-448. 

Hinze,  J.O.  1975.  T urbu I ence , second  edition.  McGraw-Hill,  790  pp. 

Kays,  W.M.  1966.  Convective  Heat  and  Mass  Transfer.  McGraw-Hill, 


i 


rr 

I 


- 115- 

Kramers,  H.  1946.  Heat  transfer  from  spheres  to  flowing  media. 

Physics,  J_2 : 61-80. 

Lighthill,  M.J.  1954.  The  response  of  laminar  skin  friction  and  heat 
transfer  to  fluctuations  in  the  stream  velocity.  Proceedings 
of  the  Royal  Society  of  London,  Series  A,  vol.  224,  pp.  1-23. 

Lueck,  R.G.  1979.  The  unsteady  flux  of  heat  and  momentum  from  a flat 
plate.  Submitted  to  Journal  of  Fluid  Mechanics. 

Morrison,  G.L.,  A.E.  Perry  and  A.E.  Samuel.  1972.  Dynamic  Calibration 

of  Inclined  and  Crossed  Hot  Wires.  Journal  of  Fluid  Mechanics, 
52:  405-474. 

Nasmyth,  P.W.  1970.  Oceanic  Turbulence.  Ph.D.  Thesis,  Institute  of 
Oceanography,  University  of  British  Columbia. 

Nowell,  Arthur  R.M.  1974.  Some  response  characteristics  of  parabolic 
hot  films  in  water.  Journal  of  Hydronautics  8_:4,  pp.  169-171. 

Schlichting,  H.  1955.  Boundary  Layer  Theory.  Pergamon  Press,  535  pp. 

Tan-atichat,  J.,  H.M.  Nagib  and  J.W.  Pluister.  1973.  On  the  interpre- 
tation of  the  output  of  hot-film  anemometers  and  a scheme  of 
dynamic  compensation  for  water  temperature  variation.  Symposium 
on  Turbulence  in  Liquids,  University  of  Missouri,  Rolla. 

White,  F.M.  1974.  V i scous  Fluid  Flow.  McGraw-Hill,  725  pp. 


This  paper  could  not  have  been  produced  without  the  co- 


operation of  the  Ocean  Mixing  Group  at  the  Institute  of  Ocean  Sciences, 
Patricia  Bav,  particularly  Dr.  A.  Gargett,  Dr.  P.  Nasmyth  and  Mr.  G. 
Chase.  Mr.  D.  Evans,  formerly  with  the  Defense  Research  Establishment, 
Pacific,  has  provided  needed  information  on  the  probes'  electronics. 
Prof.  T.R.  Osborn  deserves  appreciation  for  his  encouragements  and 
support.  This  work  has  been  supported  by  the  Office  of  Naval  Research, 
Contract  Number  00014-083-207. 


117 


118- 


Abstract 

A model  for  the  heat  flux  from  coated  spherical  and  plate- 
shaped thermistors  to  their  fluid  environment  has  been  developed 
and  agrees  favourably  with  the  available  data.  The  heat  flux  model 
is  used  to  derive  the  thermistor's  quasi-steady  sensitivity  to  speed 
and  to  temperature.  The  coating  reduces  significantly  both  the 
sensitivity  to  speed  and  to  temperature  as  well  as  the  ratio  of 
speed-to-temperature  senstivity.  The  available  values  for  the 
rates  of  temperature  and  kinetic  energy  dissipation  in  the  ocean 
indicate  that  (1)  velocitv  measurements  bv  heated  thermistors  may 
be  contaminated  by  temperature  in  the  main  thermocline  and  (2) 
temperature  measurements  by  a nearly  unheated  thermistor  may  be 
contaminated  by  velocity  in  the  mixed  laver. 
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1.  Introduction 

Thermistors  are  now  commonly  used  to  measure  temperature  micro- 
structure in  the  ocean.  Their  large  sensitivity  to  temperature,  due 
to  a high  fractional  temperature  coefficient  of  resistance  (approxi- 
mately 0.04/°C)  and  a large  resistance  (up  to  25  x 106  ohms),  make 
thermistors  attractive  for  oceanographic  use.  There  are,  however, 
several  problems  with  using  thermistors.  Three  of  these  problems  are 
(1)  the  thermistor's  signal-to-noise  ratio  which  has  been  examined 
by  Gregg,  Meagher,  Pederson,  and  Aagaard  (1978),  (2)  its  frequency 
(and  spatial)  response  which  has  been  examined  by  Lueck,  Hertzman, 
and  Osborn  (1977),  and  (3)  the  sensitivity  to  velocity  due  to  self- 
heating which  has  been  raised  by  Gregg,  et  al.  (1978). 

Inevitably,  a small  but  finite  electrical  current  must  be  passed 
through  the  thermistor  to  measure  its  electrical  resistance  and  hence 
its  temperature.  The  electric  power  dissipated  by  the  ohmic-heating 
raises  the  temperature  of  the  thermistor  above  the  ambient  fluid's 
temperature  and  thereby  subjects  the  thermistor  to  a velocity-dependent 
forced  convective  cooling. 

At  a given  resistance,  the  sensitivity  to  temperature  is  propor- 
tional to  the  current  passing  through  the  thermistor.  Hence,  the  signal 
to-noise  ratio  is  readily  increased  by  increasing  the  thermistor  current 
The  limit  to  which  the  current  can  be  increased  is  governed  by  the 
thermistor's  relative  velocity-to-temperature  sensitivity,  the  relative 
velocity-to-temperature  activity  in  the  fluid  being  probed,  and  the 
level  of  velocity  contamination  in  the  temperature  signal  that  one  is 
willing  to  accept. 
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A deliberate  large  self-heating  is  also  interesting  because  it 
turns  the  thermistor  into  an  anemometer  with  some  level  of  temperature 
contamination.  The  use  of  thermistors  as  anemometers  has  so  far  been 
limited.  Caldwell  (1968)  reported  on  an  indirectly  heated  thermistor 
to  measure  bottom  currents.  Lange  (personal  communication)  has  mounted 
both  0.015  cm  micro-bead  and  0.05  cm  flake  thermistors  on  a free-fall 
vehicle  for  microstructure  measurements.  The  thermistors  were  heated  by  a 
constant  current  source  to  approximately  50°C  above  the  ambient  temperature. 

Their  relatively  large  temperature  coefficient  and  large  resistance 
compared  to  thin  metal  film  sensors  also  make  thermistors  attractive 
for  measuring  velocity  microstructure.  In  this  paper  we  propose,  in 
section  2,  a thermal  model  for  the  heat  flux  from  thermistors.  Section 
3 compares  the  model  with  the  available  data.  The  quasi-steady  sensi- 
tivity of  the  thermistor  to  speed  and  temperature  in  the  ambient  fluid 
is  derived  in  section  4 both  for  constant  temperature  operation,  4a, 
and  constant  current  operation,  4b.  The  implications  for  thermometry 
and  anemometry  in  the  ocean  are  discussed  in  section  5. 
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2.  Steady  Response 

Following  Lueck,  et  al.  (1977)  we  will  use  a sphere  to  model 
microbead  thermistors  and  a flat  plate  for  flake  thermistors.  Both 
models  are  shown  in  Figure  1.  The  models  consist  of  an  inner  core 
of  metal-oxide,  a uniform  coating  of  thickness  A and  a boundary  layer 
of  unspecified  thickness.  Because  the  ohmic  heating  in  the  metal- 
oxide  is  distributed  over  the  volume  of  the  core  we  will  assign  a 
spatially  uniform  temperature  T to  the  metal-oxide  core.  The  temper- 
ature of  the  wetted  surface  is  0^  and  the  temperature  in  the  ambient 
fluid,  at  the  outside  limit  of  the  boundary  layer,  is  9. 

2a.  Sphere 

The  heat  flux,  Q,  from  a sphere  with  a uniform  surface  temperature 

0 is 
s 

Q = 2tt(R  + A)X2  (0s  - 9 )Nu  (1) 

where  Nu  is  the  Nusselt  number,  X2  the  fluid's  thermal  conductivity,  and 
2(R  + A)  is  the  diameter  of  the  sphere.  In  the  coating  the  temperature 
distribution,  if  spherically  symmetric,  is  governed  by 


which  is  satisfied  by 

P + A 

T (r)  = T - (T  - 0 ) —t — 

c s A 


R < r < R + A 
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Under  steady  conditions  the  heat  flux  through  the  coating  must  equal 
the  heat  flux,  Q,  into  the  fluid.  Hence,  we  use 

= Q 


3 Tc 


-A,  4tt  R2 
1 3 r 


R 


where  Aj  is  the  thermal  conductivity  of  the  coating,  to  eliminate  tr 
surface  temperature,  9 , and  derive 

T - 6 A/P  1 

Q 4f ( R + A)Aj  2it(P  t A)A2Nu 

Eauation  2 relates  the  heat  flux  from  the  model  ■thermistor  to  the 
temperature  difference  (T  - 9),  to  the  properties  of  the  coa^i.  g and 
the  properties  of  the  flow.  If  the  ratio  (T  - 0)/O  is  regarded  as  . -h>-n  1 
resistance,  in  analogy  to  Ohm's  Law,  where  T - 9 is  the  thermal  "voltage" 
and  Q the  thermal  "current",  then  the  first  term  on  the  righc  side  in 
equation  2 represents  the  thermal  resistance  of  the  coating  ani  the 
last  represents  the  thermal  resistance  of  the  boundary  layer.  The 
thermistor's  inner  core  thus  "sees"  two  thermal  resistances  between  i:: 
and  the  ambient  fluid,  namely  the  coating  and  the  boundary  layer. 

The  velocity  dependence  of  the  heat  flux  is  implicit  in  the 
Nusselt  number.  The  Nusselt  number  for  spheres  of  approximately 
uniform  surface  temperature  was  empirically  determined  bv  Kramers  ( lrug } 
as 

Nu  = 2.0  + 1.3  Pr°’15  + 0.66  Pr°'31  Re1/2  (3) 
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where 


2(R  + A) 


Re  is  the  Reynolds  number,  Pr  the  Prandtl  number,  U the  ambient  fluid's 
speed,  V the  ambient  fluid's  kinematic  viscosity,  and  K its  thermal 
diffusivity. 

The  effect  of  the  coating  on  the  thermistor's  heat  flux  is  more 
readily  determined  by  rearranging  equation  2 to  read 


2tt(P  t A)A2  (T  - 9)NU 

f A2_  ANu 
11  Xi  J 


For  an  uncoated  sphere  A = 0 and  hence  the  heat  flux  relates  to  the 
fluid's  speed  in  a fashion  after  King's  Law  because  of  the  Nusselt- 
Reynolds  number  relationship  shown  in  equation  3.  A finite  coating 
thickness,  however,  reduces  the  heat  flux's  dependence  on  speed  through 
the  term 


_ X2_  ANu 

‘ Xi  2R 


The  term  DF  depends  on  the  ratio  of  the  fluid-to-coating  conductivity, 
the  relative  coating  thickness,  A/2R,  and  the  Nusselt  number  of  the 
flow.  It  can  also  be  regarded  as  the  ratio  of  the  coating-to-boundary 
layer  resistance.  Evidently,  DF  is  an  important  design  parameter  because 
if  DF  « 1,  there  is  no  coating  effect,  and  if  DF  >>  1,  the  heat  flux  is 
independent  of  the  flow.  In  addition,  equation  4 shows  that  the  ratio 


"T' 


(T  - 0)/Q  is  a constant  at  a given  Nusselt  number  and  that  the  heat  flux 
does  not  follow  King's  Law  if  DF  is  not  small  compared  to  unity. 

The  term  DF  can  be  quite  large  for  commonly  available  thermistors.  For 
example,  0.015  cm  appears  to  be  the  lower  limit  for  the  diameter  of  the 
metal-oxide  core  (M.  Sapoff,  personal  communication)  and  at  a speed  of  1C 
cm/sec  the  Nusselt  number  is  approximately  10  in  water.  Accordingly,  a 
0.015  cm  diameter  microbead  coated  with  paralene-C  = 0.0010  watts/"* C cm) 

to  a thickness  of  A = 0.0018  cm  has  a DF  value  of  7.  The  same  thermistor 
under  the  same  flow  conditions  but  coated  with  0.0025  cm  of  glass  (Xj 
0.0105  watts/°C  cm)  has  a DF  value  of  0.9. 

2b.  Plate 

The  total  heat  flux  from  both  wetted  surfaces  of  a plate  held  at  a 

uniform  temperature  0 is  (Schlichting , 1955  , p.  269  ) 

Q = 2bX .( 0 - 0 )Nu 

2 s 

where 

1/3  1/2 

Nu  = 0.664  Pr  Re 


v 


b is  the  plate's  width  and  L its  length  in  the  flow  direction.  The  heat  flux 

through  the  coating  is 
2bLA 

Q = — r-i-  (T  - 0 ) (6) 

A s 

The  surface  temperature  is  readily  eliminated  to  yield  the  plate's  thermal 
resistance  as 

= Mk  + L_  . (7> 

0 2b  X , 2bX„Nu 
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3.  Observations 
3. a Plate 

Lange  (personal  communication)  measured  the  temperature  of  a flake 
thermistor  for  various  heating  currents.  The  flake  thermistor  had 
the  nominal  dimensions  of  0.051  by  0.051  by  0.0025  cm  and  was  coated 
with  paralene-C  to  a nominal  thickness  of  0. 0018  cm.  The  measurements 
were  made  by  towing  the  thermistor,  with  its  face  parallel  to  the 
flow,  down  the  length  of  a 30  by  2.4  by  0.24  metre  tow  tank  at  speeds 
of  10,  20,  30,  40,  50  and  60  cm/sec.  The  accuracy  of  the  speed  measure- 
ments is  15%  or  better.  We  used  Lange's  data  to  calculate  the  heat 
flux  and  the  thermistor's  (metal-oxide's)  temperature  for  the  various 
heating  currents  and  flow  speeds.  The  results  are  shown  in  Figure  2. 

The  slopes  of  the  temperature  vs  heat  flux  plots  for  a given  speed  equal 

the  thermal  resistance  (T  - 9)/Q  (because  T = 9 when  Q = 0).  The 

-1/2  . 

thermistor's  thermal  resistance  is  plotted  against  U m Figure  3. 

The  observed  relationship,  by  a least  square  fit,  is 


0.437  ± 0.001 


0,538  t 0.005 

^1/2  mw 


(8) 


Using  the  appropriate  parameter  values  listed  in  Table  1 yields, 
according  to  equation  7, 


T - 9 


Q 


0.35 


0.59 

+ ^1/2  mw 


The  observed  and  the  predicted  thermal  resistances  agree  favorably 
considering  that  the  flake's  dimensions  are  known  only  by  the  nominal 
values  specified  by  the  manufacturer. 
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Figure  2.  Heat  flux  versus  temperature  for  various  speeds  for 
Lange's  0.051  x 0.051  x 0.0025  cm  paralene-c  coated 
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3b.  Sphere 

Using  the  equipment  specified  above,  Lange  (personal  communication) 
also  tested  an  ellipsoidal  microbead  thermistor  with  a nominal  semi- 
minor axis  diameter  of  0.013  cm  coated  with  paralene-C  to  a thickness 
of  0.0018  cm.  The  temperature  vs  heat  flux  curves  for  the  various 
tow  speeds  are  shown  in  Figure  4.  We  have  calculated  the  microbead's 
Nusselt  number  using  equation  3 and  the  parameter  values  in  Table  1. 

The  microbead's  thermal  resistance  is  plotted  against  the  inverse  of 
its  Nusselt  number  (in  Figure  5).  The  least  square  fitted  line  in 
Figure  5 is 


0.59  + 


2.51 
Nu  mw 


In  comparison,  the  relationship  predicted  by  equation  2 is 


2.85 


3.35 
Nu  mw 


The  agreement  between  the  predicted  and  the  observed  coating  resistance 
(2.85  vs  0.59  °C/mw)  is  poor  and  the  agreement  between  the  predicted 
and  the  observed  Nusselt  number  coefficient  (3.35  vs  2.51)  is  only 
fair.  However,  in  the  absence  of  more  precise  knowledge  about  the 
thermistor's  diameter  and  its  coating  thickness  it  is  not  clear  if 
the  thermal  model  is  in  error  or  if  the  parameters  used  are  erroneous. 

For  further  tests  we  ordered  a set  of  20  microbead  thermistors 
from  one  production  batch.  The  thermistors'  semi-minor  axis  are, 
nominally,  0.015  cm.  Half  (10)  of  the  thermistors  were  coated  with 


glass  to  a nominal  thickness  of  0.00254  cm.  All  20  microbeads  were 
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photographed  through  a microscope  to  determine  their  outside  dimensions. 
The  difference  between  the  dimensions  of  the  coated  and  the  uncoated 
microbeads  was  used  to  derive  a statistical  estimate  for  the  average 
coating  thickness.  The  derived  coating  thickness  is  0.0028  +_  0.007 
cm  which  is  only  about  10%  larger  than  the  nominal  value.  A back 
lighted  microphotograph  of  the  two  coated  thermistors  tested,  namely 
#9  and  #10,  is  shown  in  Figure  6. 

The  tank  used  to  test  the  thermistors  is  an  annulus  of  width 
4 cm,  depth  4 cm  and  outer  radius  14  cm.  The  annulus  was  placed  on 
a rotating  table.  The  speed  of  the  fluid  past  the  stationary  thermistor 
was  determined  from  measurements  of  the  rotation  rate  and  the  radial 
position  of  the  thermistor  under  the  assumption  of  solid  body  rotation 
in  the  fluid.  The  walls  of  the  annulus  were  roughened  to  reduce  fluid 
motions  relative  to  the  rotating  table. 

Samples  of  the  temperature  vs  heat  flux  relationship  for  the 
thermistors  are  shown  in  Figure  7.  The  thermal  resistance  is  plotted 
against  the  inverse  of  the  Nusselt  number  in  Figure  8.  It  should  be 
noted  that  each  point  on  these  plots  represents  10  observations  of 
the  thermistor’s  temperature  and  heat  tiux.  The  Nusselt  number  used 
is  based  on  equation  3 and  the  parameters  in  Table  1.  One  data  point 
in  Figure  8(b)  is  at  variance  with  the  straight  line  fitted  through 
the  other  four  points . This  data  point  comes  from  the  first  flow 
speed  tested.  The  point's  deviation  cannot  be  explained  as  a 
"blunder"  in  either  the  speed  measurement  (33.6  cm/sec)  or  the  heat 
flux  vs  temperature  regression.  The  consistency  of  the  zero  heat 
flux  intercept  (Fig.  9b)  for  various  speeds  also  rules  out  a sporadic 
change  in  the  electrical  resistance  temperature  relation  for  the 
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Figure  6 


Back  lighted  micrograph  of  the  authors’  glass  coated 
microbeads  and  a sketch  of  the  thermistors  mounting 
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Figure  7(a)  Heat  flux  versus  temperature  for  various  flow  speeds  for 
thermistor  number  9. 
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thermistor.  The  point's  variance  is  most  likely  due  to  a surface 
contamination,  such  as  lint,  which  was  too  small  to  be  noticed  during 
the  measurements  and  dislodged  prior  to  the  readings  at  other  speeds. 
Breakage  of  the  probe  precludes  a re-examination  of  the  thermistor's 
heat  flux  at  that  speed.  Because  the  other  thermistors  tested  do  not 
have  such  a deviant  data  point  we  have  not  used  it  to  fit  the  line 
shown  in  Figure  8(b).  We  have,  however,  retained  the  data  point 
because  we  do  not  have  an  explicit  reason  for  removing  it. 

The  least  sauare  fitted  thermal  resistances  are 


T - 9 


Q 


0.527  ± 0.001 


2.09  t 0.01 

Nu  mw 


for  microbead  #10  and 


(9) 


0.422  ± 0.001 


1.92  ± 0.02  ^C_ 

Nu  mw 


for  microbead  #9.  In  comparison,  the  predicted  thermal  resistance, 
using  equation  2,  is 


0.21  + 


2.4 

Nu 


°C 

mw 


for  both  thermistors.  The  observed  Nusselt  number  coefficients  for  the 
two  thermistors  are  in  reasonable  agreement  with  the  predicted  values. 
The  observed  thermal  resistance  of  the  coating  is,  however,  2 to  2.5 
times  larger  than  predicted  and  consequently  warrants  some  consideration 
The  model  for  the  coating's  temperature  is  based  upon  spherical 


symmetry.  The  fact  that  the  thermistors  are  ellipsoidal  rather  than 
spherical  cannot  alone  account  for  the  observed  discrepancy.  However 
it  is  known  that  the  heat  flux  from  bluff  bodies  is  preferentially 
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biased  toward  their  front  surface  where  the  boundary  layer  is  thinnest. 

Thus  the  surface  through  which  a large  fraction  of  the  heat  passes  is 
much  smaller  than  the  total  surface  area  of  the  sphere.  Perhaps  the 
front  hemisphere,  or  even  less,  represents  the  area  that  transmits  the 
heat  from  the  core  to  the  ambient  fluid.  An  examination  of  equation  2 
then  indicates  that  the  term  representing  the  coating's  resistance 

A 

4ttR(R  + A )X] 

should  really  be  larger  to  account  for  this  effective  surface  area 
reduction.  Until  the  heat  flow  process  is  better  understood,  multi- 
plying the  above  coating  resistance  by  a factor  of  about  2 or  2.5 
might  maintain  the  usefulness  of  equation  2. 

The  Nusselt-Prandtl-Reynolds  number  relationship  for  spheres,  on  the 
other  hand,  is  an  empirical  observation  independent  of  the  local  heat 
flux  distribution  over  the  sphere's  surface.  The  difficulty  with  using 
Kramers'  (1946)  Nusselt  number  formulation  (equation  3)  is  that  it 
applies  only  to  spheres  with  a nearly  uniform  surface  temperature. 

Because  the  heat  flux  from  a sphere  is  not  uniform  over  its  surface, 
the  presence  of  the  coating  leads  to  a non-uniform  temperature  over 
the  wetted  surface.  Even  for  simple  geometries  such  as  a plate,  the 
effect  of  a non-uniformity  of  the  surface  temperature  on  the  heat 
flux  is  difficult  to  predict.  Calculations  by  White  (1974)  suggest  that 
the  coating  should  make  the  heat  flux  more  uniform  over  the  wetted  surface. 
The  fact  that  the  observed  Nusselt  number  coefficients  differ  from  the 
predicted  coefficients  by  only  13  and  20%  indicates  that  Kramers' 
relationship  has  some  validity  for  the  thermistors  tested. 
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There  is,  however,  a more  general  problem  with  predicting  the 
response  of  thermistors.  There  is  evidence  for  non-uniformity  between 
nominally  identical  thermistors.  This  is  particularly  evident  in  the 
coating  thickness  which  may  vary  considerably  over  the  thermistor's 
surface  and  between  thermistors,  as  noted  by  Gregg  et  al.  (1978)  and  by 
us.  Although  we  were  not  able  to  observe  directly  the  coating  thicknesses 
of  the  particular  thermistors  reported  here,  the  difference  between  the 
observed  response  of  thermistors  #9  and  #10,  whose  outside  dimensions 
differed  by  less  than  2%,  suggest  that  their  coatings  may  be  quite 
different  in  thickness  and  in  thickness  distribution.  Because  of  this 
non-uniformity  among  nominally  identical  thermistors  it  is  not  clear 
that  improvements  of  the  thermal  model,  although  intellectually 
interesting,  can  lead  to  corresponding  improvements  in  response 
predictability. 
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4.  Quasi-Steady  Response 

The  response  is  considered  to  be  quasi-steady  when  the  heat  flux  at 
any  instant  is  that  appropriate  to  steady  conditions  at  the  instantaneous 
value  of  the  fluid's  speed  and  temperature;  i.e.  the  response  in  the 
limit  of  zero  frequency. 

4a.  Constant  Temperature  T 
(i)  Sphere 

The  variation  of  equation  2 with  respect  to  speed  and  temperature 

gives 

T - 9 dQ  d§  0.33  Pr°‘  ilRe1/2  dU 

Q Q Q 2tt(R  t A)X2Nu2  U ( 

The  sensitivity  of  the  heat  flux  to  temperature  is 


dQ 

d0 


dU=0 


-Q 

T - 6 


(11) 


which  is  a constant  at  a given  speed  and  readily  determined  from  the 
temperature  vs  heat  flux  relationship.  For  thermistor  #10,  for  example, 
the  sensitivity  to  temperature  is  1.39  mw/°C  at  15.4  cm/sec.  The 
sensitivity  of  the  heat  flux  to  speed  is 


dQ 

f Q 1 

0.33  Pr‘  '1Re1/2 

T - 9 

dU 

d6=0 

-3 

1 

O 

[2tt ( R + A)X2Nu2  J 

U 

which  at  15.4  cm/sec,  using  only  the  dimensional  information  in  Table  1, 
equals  0.0111  mw/cm/sec  for  every  centigrade  degree  of  overheat  (T  - 0). 
Choosing,  say,  T - 9 = 50°C  results  in  a sensitivity  of  0.553  mw/cm/sec. 
The  relative  speed-to-temperature  sensitivity  then  is 
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(dQ/dU )d0=O 

f Q 1 

0.33  Pr,31ReI/2 

T - © 

CdQ/d0)dU=O 

t - e 

l2tt(R  + A)X2Nu2 

U 

and  equals  0. 40°C/cm/sec  for  thermistor  #10  under  the  conditions 
specified,  namely  U = 15.4  cm/sec  and  T - 0 = 50°C.  The  implications 
are  that  a quasi-steadv  fluctuation  of  1 cm/sec  is  indistinguishable 
from  a fluctuation  of  0.4  °C.  Alternate  expressions  for  the  sensi- 
tivities that  clarify  the  effect  of  the  coating  are 


dQ 

dU 


d6=0 


1 1 1/9 

2 (R  + A)^2 [0.33  Pr’  Ra  1 T - 6 
(1  + DF)2  U 


(14) 


and 


dQ  _ 2tt(R  + A) A?Nu 

d0  dU*0  " 1 + DF 


(dQ/dU )dQ=0  _ 0.33  Pr'jlRe1/2  T - 9 
(dQ/d© )dU=0  ' (1  + DF)Nu  U 


where  DF  represents  the  coating-to-boundary  layer  resistance  ratio. 

For  thermistor  #10  the  observed  DF  value  at  15.4  cm/sec  is  2.2,  which 
implies  that  the  speed  sensitivity  has  been  reduced  by  a factor  of 
(1  + 2.2)2  = 10  from  its  value  in  the  absence  of  a coating.  The 
temperature  and  the  relative  speed-to-temperature  sensitivity,  using 
equations  15  and  16,  is  reduced  by  a factor  of  (1  + 2.2)  = 3.2  by  the 
coating.  Clearly  the  coating  requires  careful  consideration  when  one 
chooses  a thermistor  for  heated  anemometry  and  heated  thermometry. 

We  have  used  the  observed  thermal  resistance  of  thermistor 
number  10  (eq.  9)  to  predict  the  speed  sensitivity  of  that  thermistor 
for  various  flow  speeds.  The  speed  sensitivity  per  degree  of  overheat 
is,  using  equation  9 and  the  parameters  in  Table  1, 
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7 


Lange  Authors 


Plate 

Microbead 

Microbeads 

# 9 

# 10 

R 

(cm) 

— 

0.0064 

0.0087 

0.0087 

L 

(cm) 

0.051 

— 

— 

— * - 

b 

(cm) 

0.051 

— 

— 

— 

A 

(cm) 

0.0018 

0.0018 

0.0028 

0.0028 

(watts/cm°C) 

0.0010 

0.0010 

0.010 

0.010 

*2 

(watts/cm°C) 

0.0059 

0.0059 

0.0059 

0.0059 

V 

(cm2 /sec) 

0.0123 

0.0120 

0.0097 

0.010 

K 

(cm2/sec) 

0.00139 

0.00139 

0.0014 

0.0014 

Pr 

8.9 

8.63 

7.1 

6.8 

Table  Cl 

The  parameter 

value  associated  with  the 

thermistors 

examined. 

The 

dimensional  data  for  Lange's 

thermistors 

are  nominal 

values . 

The 

conductivity  of 

paralene-C  (X 

j for  Lange* 

's  thermistors)  is 

after  Alpaugh  and  Morrow  (1974). 

The  viscosity  (v),  the 

diffusivity 

(k),  and  the  Prandtl  number  (Pr)  are  after  Batchelor  (1970,  p.  597). 
The  conductivity  of  glass  for  the  authors'  thermistors)  is 
after  Lueck,  et  al.  (1977). 
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1 dQ 

T - 9 dU 


0.230/U' 


in  units  of  °C  per  cm/sec  per  degree  cent i grade,  which  is  plotted  in 

Figure  9.  Both  the  speed  sensitivity  and  the  ratio  of  speed-to- 
temperature  sensitivity  decrease  rapidly  with  increasing  flow  speed. 

Consequently  the  usefulness  of  heated  thermistors  for  anemometry  may 

be  restricted  to  low  flow  speeds. 


(ii)  Plate 


The  variation  of  equation  7 gives  the  sensitivity  of  the  plate's 


heat  flux  as 


0.664  bA,  Pr1/3Re1/2  T - 9 


(1  + y2 

A j Lj 


Q 

T - 9 


1/3  1/2 

-2bA20. 664  Pr  Re 

(1  + j]  r Su) 


(dQ/dU)d9=0 

(dQ/d9)dU=0 


1 A2A.. 

1 + 17  t "u 


T - 9 
2U 


Using  the  observed  thermal  resistance  shown  in  equation  (8)  we  calculate 
the  relative  sensitivity  as 


(dQ/dU)d9=0  _ 1.00  T-0 

(dQ/d9)dU=0  1 + 0.825U1/2  2U 


which,  for  comparison  with  the  glass-coated  microbead  #10  equals 
0.38  °C/cm/sec  at  T - 9 = 50  °C  and  U = 15.4  cm/sec.  The  term  DF 

I/O 

equals  0.825  U for  the  plate  tested  by  Lange. 


AMBIENT  FLUID  SPEED  U cm/sec 


of  overheat  for  fhe  authors'  g(ass  coated  thermistor 
number  10  and  Lange's  paralene-c  coated  flake  thermistor. 
The  curves  are  based  on  the  observed  thermal  resistances 
given  by  equations  8 and  9. 


Do  ep/op  e-i 

oas/uuD/Oo  OP/OP  I 
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The  relative  speed-to-temperature  sensitivity  of  Lange's  flake 
thermistor  is  plotted  in  Figure  10.  The  speed  sensitivity  per  degree 


of  overheat  is  given  by 


1 dQ 
T - 0 dU 


1.86 

(1  + 0. 812*701  2U 


which  is  plotted  in  Figure  9. 

The  assymptotic  value  of  the  speed  sensitivity  of  a heated  sensor 

-1/2  -3/2 

at  high  speeds  is  proportional  to  U without  a coating  and  U 

with  a coating.  The  assymptotic  value  of  the  ratio  of  the  speed-to- 

— 3/2 

temperature  sensitivity  is  U-1  without  a coating  and  U with  a 
coating . 


4b . Constant  Current  I 
(i)  Sphere 

In  the  constant  current  mode  the  temperature  of  the  thermistor's 
core  is  variable  and  is,  in  practice,  the  sensed  quantity.  The  variation 
of  the  heat  flux  dQ  is  related  to  the  variation  of  the  metal-oxide’s 
temperature  by 


dQ  = IdE  = I dRT  = QadT 


where  I is  the  constant  current  passing  through  the  thermistor,  dE  is  the 
fluctuation  of  the  thermistor's  voltage,  dR^  the  fluctuation  in  its 
electrical  resistance,  dT  the  fluctuation  in  its  temperature,  and  a its 
fractional  temnerature  coefficient  of  electrical  resistance  defined  by 


1 dRT  ~ 


R?  dT 


= -0 . 04/°C 


The  variation  of  equation  2 is 


dT  d9  _ T - 9 dO 
Q ' Q Q 0 


O I I/O 

0.33  Pr* ' Re  dU 

2tt(P  + A)X,Nu2  U 


which,  when  combined  with  equation  20,  yields  for  the  sphere, 

31  1/2 

r,  ../*  ^-.dT  d9  _ 0.33  Pr*  Re  7 dU 
[!  _ a(T  - 9)]—  - — - - - A)Mu2  — 


The  constant  current  temperature  sensitivity  is 


1 - oc(T  - 9) 


According  to  (22)  the  temperature  sensitivity  is  unity  at  zero  overheat 
which  is  desired  if  the  thermistor  is  to  be  used  as  a temperature 
sensor  with  infinitesimal  internal  heating.  The  temperature  sensi- 
tivity decreases  with  increasing  internal  heating  (noting  that  a < 0) 
and  vanishes  at  very  high  levels  of  internal  heating.  It  is  very 
interesting  to  observe  that  sensors  with  positive  temperature  coeffi- 
cients (a  > 0)  increase  their  temperature  sensitivity  with  increasing 
internal  heating  in  the  constant  current  mode.  The  condition  T - 0 = a-1 
corresponds  to  an  instability  for  positive  temperature  coefficient 
sensors.  In  contrast,  the  temperature  sensitivity  in  the  constant 
temperature  mode  is  independent  of  the  overheat  T - 9. 

The  speed  sensitivity  of  the  spherical  thermistor’s  heat  flux  is 


Q 0.33  Pr,31Re1/2 

T - 9j  [1  - a(T  - 9 ) ] 2tt ( P + A)A,NuJ 


T - 9 
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At  zero  overheat  (T  = 9)  the  thermistor  has  no  speed  sensitivity. 

The  speed  sensitivity  increases  with  increasing  overheat  to  the  limit 


dT  1 0.33  Pr'31Re1/2  Q 1 

dU  ..  a 2tt(R  + A) A, Nu2  T - 9 U 

d9=o  2 


limit  Q -*■  00 


which  is  independent  of  the  overheat  because  Q/(T  - 9)  is  a constant 
at  a given  Nusselt  number.  The  ratio  of  the  quasi-steady  speed-to- 
temperature  sensitivity  is 


(dT/dU)d9=0 
( dT/d9 )dU=0 


Q 0.33  Pr'31Re1/2  T - 9 
T - 9j  2tt(R  + A)A2Nu2  U 


which  is  identical  to  the  relative  sensitivity  in  the  constant 
temperature  heating  mode.  The  advantages,  if  any,  of  operating  the 
thermistor  in  the  constant  temperature  mode  over  operating  it  in  the 
constant  current  mode  are  reflected  in  the  unsteady  response  which 
has,  so  far,  not  been  measured. 
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5 . Discussion 

Some  knowledge  of  the  relative  activity  of  temperature  and  velocity  in 
the  fluid  being  probed  (such  as  the  ocean)  is  necessary  to  determine  the 
suitability  of  heated  thermistors  for  turbulence  measurements  and  to  assess 
the  velocity  contamination  of  (nearly)  unheated  thermistors  used  for  thermo- 
metry. We  will  use  as  a basis  for  discussion  some  of  the  available  observa- 
tions of  the  kinetic  energy  dissipation  rate  e and  the  temperature  dissipation 


rate  x where,  assuming  isotropy 


e = 15v 


' ^ ] 
. 9z  J 


(24a)  , x = 


(24b) 


3w'/3z  is  the  vertical  gradient  of  the  vertical  component  of  the  velocity 
fluctuations  and  30' /3z  is  the  vertical  gradient  of  the  temperature  fluctua- 
tions. Gregg  (1976)  reported  temperature  dissipation  rates  of  4.3  x 10~8  to 
1.1  x 10" *0  °C2/sec  averaged  over  several  tens  of  metres  in  the  mixed  layer 
and  values  of  3.6  x 10-8  to  2.1  x 10-8  °C2/sec  averaged  over  several  tens  of 
metres  in  the  thermocline  of  the  central  North  Pacific.  Gregg  (1977)  also 
observed  temperature  dissipation  values  of  2.2  x 10-12  to  6 x 10"11  °C2/sec 
for  depths  below  1000  metres  in  several  widely  spaced  locations  in  the  Pacific. 
Gargett  (1978)  measured  dissipation  rates  as  large  as  9.8  x 10-7  °C2/sec  in 
an  upper  ocean  frontal  regime. 

The  rate  of  kinetic  energy  dissipation  near  the  Azores  has  been  measured 
by  Osborn  (1978).  Averaged  over  5 metres  of  depth  the  values  were  typically 
10” 14  cm2 /sec 8 but  ranged  as  high  as  10"2  cm2 /sec8.  Crawford  and  Osborn  (1978) 
reported  averaged  kinetic  energy  dissipation  rates  of  up  to  7 x 10  8 , 

4 x 10" 5 and  2 x 10” 8 cm2 /sec8  for  the  regions  above,  in  and  below  the  core 
of  the  Atlantic  Equatorial  Undercurrent. 
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If  we  denote  the  ratio  of  a probe's  speed-to-temperature  sensi- 
tivity by  B,  then  the  velocity  contamination  of  a temperature  dissi- 
pation rate  measurement  is,  using  24a  and  b. 


X 


cont 


(25) 


Correspondingly,  the  temperature  contamination  of  a kinetic  energy 
dissipation  rate  measurement  is 


5Pr  1 

econt  2 6^"  X 


(26) 


Using  Gargett's  (1978)  value  of  9.8  x 10-7  °C2/sec  as  an  upper 
limit  for  "X  in  the  thermocline,  and  using  the  observed  relative  sensi- 
tivity of  0.29  °C/cm/sec  for  thermistor  #10  (Fig.  10)  at  T - 9 = 50°C, 
U = 15.4  cm,  and  Pr  = 8 gives  a contamination  of 

£ = 2.3  x 10-1*  cm2/sec3 

cont 

which  is  comparable  to  the  values  of  e reported  by  Osborn  (1978)  in 
the  seasonal  thermocline  adjacent  to  the  island  of  Santa  Maria  in  the 
Azores.  In  the  mixed  layer,  using  Gregg's  (1976)  upper  value  of 
4.3  x 10-9  °C2/sec , 


£ = 1.0  x 10-6  cm2/sec3 

cont 


which  may  be  small  compared  to  the  kinetic  energy  dissipation  in  the 
mixed  layer.  At  large  depth  the  contamination  is,  using  x = 6 x 10“ 11 
°C2/sec,  1.4  X|(3-8  cm2/sec3.  However,  there  are  no  values  of  e avail- 


If  therm!  -tor  #1  w-  : 


t-  r . • 3ture  :•  ■ u uj  • m-  ; 


then  the  level  of  velocity  contaminati  n in  th<  - • : ured  temperatur< 

dissipation  rate  1;  from  equation  25 

y . =1.7  x 10”6  (T  - 9)2e 
cont 

where  we  have  used,  from  fir ure  10, 

MU,  10-3  <T  . ,) 

at  U = 15.4  cm/sec.  Taking,  for  example,  an  upper  value  of  10”2 
cm2  /sec3  for  e in  the  mixed  layer,  and  x - 10”12  °C2/sec  as  the 
lowest  dissipation  rate  to  be  resolved  implies  a 1%  contamination 
at  an  overheat  (T  - 6)  of  only  7.7  x 10”  3 °C.  Using  Figure  7b  we 
calculate  that  this  temperature  rise  results  from  only 


0 = 1.39  ^x  7.7  x 10‘3  °C  = l.l  x io“s  watts 


of  internal  heating  in  the  thermistor.  Gregg  et  al.  (1978)  reports 
internal  heating  levels  of  4 to  13  x 10” p watts  for  his  thermistors 
which,  however,  have  a different  mounting  and  are  2.5  times  larger 
than  the  microbeads  discussed  here.  Several  people  have  suggested 
warmed  thermistors  for  high  frequency  thermometry,  however,  because 
even  very  small  amounts  of  self  heating  may  lead  to  velocity  contamination 
it  appears  that  grossly  and  moderately  heated  thermistors  cannot  be  very 
useful  for  thermometry  in  the  ocean. 

We  have  used,  in  the  above  discussion,  temperature  dissipation 
rates  that  were  averaged  over  several  tens  of  metres  or  more  and 
kinetic  energy  dissipation  rates  that  were  averaged  over  5 metres 
or  more.  On  a more  local  scale,  however,  the  dissipation  rates  are 
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larger  than  cited,  and  consequently  the  amount  of  signal  contamina- 
tion may  be  more  severe  than  the  values  given  above  because  tempera- 
ture and  velocity  activity  are  not  necessarily  correlated.  We  have 
also  ignored  the  effect  of  the  heated  thermistor's  frequency  response. 
Lueck  (1979)  has  shown  that  the  frequency  response  bandwidth  of  a 
heated  plate  is  larger  for  speed  than  for  temperature.  However, 
in  the  ocean  the  wavenumber  range  of  the  temperature  gradient  spectrum 
is  typically  larger  than  the  range  of  the  velocity  gradient  spectrum. 


/ 
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6 . Conclusions 

A model  for  the  heat  flux  from  internallv  heated  spherical  and 
plate-shaped  thermistors  to  their  fluid  environment  has  been  derived. 

The  model  incorporates  the  effects  of  a thermistor's  coating  and  the 
effects  of  forced  convection  in  its  boundary  layer  through  the  use  of 
an  established  Nusselt  number.  The  model  agrees  favourably  with  the 
heat  flux  from  a paralene-C  coated  flake  thermistor  as  observed  by 
Lange.  Agreement  between  the  model  and  our  observations  for  the  heat 
flux  from  two  glass  coated  microbead  thermistors  is  also  favourable 
except  that  the  effect  of  the  coating  is  approximately  two  times  larger 
than  predicted.  Preferential  heat  conduction  through  the  frontal 
surface  of  the  thermistor  could  account  for  this  difference.  The 
model  agrees  well  with  the  observations  if  the  coating's  surface  area 
is  reduced  to  an  effective  surface  area. 

The  heat  flux  from  heated  sensors,  in  general,  depends  on  both  the 
temperature  and  the  speed  of  the  ambient  fluid.  The  simultaneous 
sensitivity  of  thermistors  to  speed  and  temperature  has  been  explicitly 
derived.  The  coating  reduces  both  the  thermistor's  sensitivity  to  speed 
and  to  temperature  from  its  value  in  the  absence  of  a coating.  The 
reduction  in  sensitivity  is  greater  for  speed  than  for  temperature,  and 
hence,  the  relative  speed-to-temperature  sensitivity  is  also  reduced 
by  the  coating. 

Some  implication  of  the  results  for  anemometry  and  thermometry 
in  the  ocean  have  been  discussed.  The  level  of  contamination  of  a 


\ 


measurement,  due  to  the  thermistor’s  simultaneous  sensitivity  to  speed 
and  temperature,  depends  on  the  amount  of  internal  heating  and  on  the 
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relative  intensity  of  velocity  and  temperature  turbulence  in  the  fluid 
being  probed.  Consequently  the  contamination  of  a measurement  depends 
strongly  on  one's  particular  circumstances.  Using  observed  values 
for  the  rate  of  temperature  and  kinetic  energy  dissipation  in  the 
ocean  we  conclude  that  it  is  possible  that  (1)  velocity  measurements 
made  by  heated  thermistors  may  be  contaminated  by  temperature  in 
the  main  thermocline,  that  (2)  temperature  measurements  made  by 
nearly  unheated  thermistors  may  be  contaminated  by  velocity  in  the 
mixed  layer  and  in  other  regions  of  high  velocity  but  low  temperature 
activity,  and  that  (3)  moderately  or  grossly  heated  thermistors  are  not 
useful  for  thermometry. 
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